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1.  INTRODUCTION  AND  SUMMARY 

The  ^vent  of  unstructured  grid  technology  has  greatly  reduced  grid  generation  time  for 
complex  configurations  and  has  the  potential  to  increase  solution  quality  through  techniques  like 
1  r^ently  these  benefits  have  been  limited  to  inviscid  applications  due  to 

e  difficulty  of  generating  unstructured  tetrahedral  grids  appropriate  for  viscous  computations, 
^ese  problems  require  grid  clustering  normal  to  the  surface  to  resolve  the  boundary  layer. 

tandard  unstructured  tetrahedral  grid  generation  methods  attempt  to  maintain  a  unitary  aspect 
ratio  causmg  the  grids  to  become  excessively  large  in  the  viscous  layer.  Recently,  a  wide  variety 
of  approaches  have  arisen  in  an  attempt  to  overcome  these  difficulties.  This  variety  of 
approaches  makes  it  difficult  to  select  a  direction  for  future  development  of  unstructured  grid 
rnethods.  This  report  documents  a  systematic  assessment  of  the  relative  strengths  and  weaknesses 
of  three  new  approaches  to  unstructured  grid  viscous  flow  field  simulation.  This  information  will 
help  guide  future  development  efforts  and  will  help  the  computational  fluid  dynamics  (CFD)  user 
select  an  appropriate  unstructured  grid  approach  for  a  given  application  based  on  accuracy  and 
efficiency  requirements. 

The  methods  evaluated  in  this  study  are  the  VGRID/USM3D,  MACGS/NASTD,  and 
TETMESH  /Cobalt  grid  generation/flow  solver  software  packages  developed  by  NASA  Langley 
Research  Center,  McDonnell  Douglas  Corporation  (MDC),  and  Air  Force  WL/FIMC  respectively. 
Each  of  these  methods  use  a  different  unstructured  grid  cell  topology  to  model  the  viscous  layer. 
An  illustration  of  the  unstructured  grid  types  used  with  each  approach  is  shown  in  Figure  1 .  The 
VGRID/USM3D  software  package  uses  tetrahedral  elements  throughout  the  domain  with 
amsotropic  stretching  near  the  surface.  This  approach  will  be  referred  to  as  the  tetrahedral 
approach.  The  zonal  hybrid  approach  used  by  the  MACGS/NASTD  software  package  combines 
a  structured  grid  consisting  of  hexahedral  elements  at  the  surface  with  tetrahedral  elements 
defined  away  firom  the  body.  The  TETMESH/Cobalt  software  package  uses  soni-structured 
layers  of  ^d  cells  in  the  viscous  region  where  the  layers  are  normal  extrusions  of  the  triangular 
surface  grid  resulting  in  prismatic  elements.  These  elements  merge  with  a  tetrahedral  grid 
defined  outside  the  viscous  region.  This  approach  will  be  referred  to  as  the  mixed  cell  approach. 

Each  of  these  modeling  approaches  have  a  unique  set  of  advantages  and  disadvantages. 
The  zonal  hybrid  approach,  while  being  straight  forward  to  implement,  does  not  offer  the 
geometric  flexibility  of  the  tetrahedral  and  mixed  cell  approaches.  This  approach  requires  a 
structured  grid  to  be  defined  near  the  geometry  surface  which  typically  takes  a  long  time  to 
generate.  The  triangular  surface  grids  used  in  the  tetrahedral  and  mixed  cell  methods  are 
generally  easier  to  generate  but  provide  poor  resolution  of  highly  curved  geometries.  Voliune 
grid  generation  is  much  easier  for  the  tetrahedral  approach  than  for  the  zonal  hybrid  or  mixed 
cell  methods  which  require  separate  grid  generation  techniques  in  the  viscous  and  farfield  regions. 
However,  grid  requirements  for  anisotropic  tetrahedral  grids  are  not  well  known  since  relatively 
few  tetrahedral  viscous  computations  have  been  performed. 
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Figure  1 .  Three  Flat  Plate  Unstructured  Grid  Topologies 

The  choice  of  grid  cell  shape  can  have  a  large  influence  on  flow  solution  efficiency.  For 
instance,  solution  computation  on  a  tetrahedral  grid  will  result  in  at  least  four  times  more  work 
per  node  than  a  hexahedral  grid.  However,  tetrahedral  grids  may  require  fewer  nodes  to  obtain 
the  same  level  of  accuracy  due  to  the  increased  directional  information  available  at  each  node. 
Solutions  computed  on  zonal  hybrid  grids  can  take  advantage  of  the  structured  data  storage  to 
provide  very  high  efficiencies,  but  generally  contain  more  grid  cells  than  tetrahedral  grids. 

In  order  to  identify  the  best  and  worst  attributes  of  each  method,  a  comprehensive  study 
of  the  strengths  and  weaknesses  of  these  approaches  has  been  performed.  This  study  includes 
an  assessment  of  the  relative  acciiracy,  efficiency,  robusmess  and  ease  of  use  of  each  method. 

An  essential  element  in  performing  an  objective  assessment  is  the  determination  of  the 
grid  resolution  required  by  each  method.  A  grid  that  is  too  coarse  will  adversely  affect  solution 
accmacy  while  a  grid  that  is  too  fine  will  reflect  poorly  on  the  method’s  efficiency.  In  this 
study,  grid  resolution  requirements  were  established  by  performing  a  grid  refinement  study  on 
laminar  and  turbulent  flat  plates  over  a  range  of  Reynolds  numbers.  For  each  unstructured  grid 
method  a  study  was  performed  to  determine  the  spacing  of  the  first  point  off  the  wall  and  the 
nximber  of  points  across  the  boundary  layer  required  to  obtain  an  accurate  solution. 

Once  the  grid  resolution  requirements  were  known,  the  viscous  prediction  capabilities  of 
each  method  were  evaluated  by  comparing  solutions  for  sample  test  cases.  These  cases  were 
selected  to  cover  a  wide  range  of  flow  field  physics.  Solutions  about  the  ONERA  M6  wing  were 
generated  at  two  run  conditions.  At  the  first  run  condition  of  Mach  0.7,  angle  of  attack  of  0.0 
degrees,  and  Reynolds  number  of  11.7  million  the  flowfield  is  subsonic  and  remains  attached 
across  the  entire  wing  span.  At  the  second  run  condition  of  Mach  0.85,  angle  of  attack  of  5.0 
degrees  and  Reynolds  number  of  11.7  million,  the  flow  field  contains  several  complex  viscous 
interactions  including  shock-induced  separation  over  the  outer  one  third  of  the  wing. 

For  each  method,  comparisons  of  solution  accuracy,  computational  requirements,  and  grid 
generation  time  were  made.  The  convergence  properties  are  compared  and  an  assessment  of  the 
ease  of  use  of  the  grid  generation  and  flow  solver  software  is  made. 
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2.  SOFTWARE  INSTALLATION 

In  this  section  a  short  description  of  each  of  the  three  software  packages  is  given.  The 
steps  that  were  employed  in  installing  this  software  on  various  computer  systems  are  also 
documented. 

The  VGRID  unstructured  grid  generation  system  was  developed  by  ViGYAN,  Inc.  under 
contract  from  NASA  Langley  Research  Center.  A  complete  description  of  the  method  is 
provided  in  Reference  (1).  VGRID  is  based  on  the  Advancing-Front  Method  for  the  generation 
of  Euler  tetrahedral  grids  and  the  Advancing-Layers  Method  for  producing  thin-layered 
tetrahedral  grids  suitable  for  computation  of  viscous  fluid  flows.  A  multi-directional  grid 
stretching  substantially  decreases  the  final  grid  size.  A  graphic  user  interface  tool,  Grid  Tool, 
is  available  for  converting  given  Computer  Aided  Design  (CAD)  geometries  into  Non-Uniform 
Rational  B-Splines  (NURBS)  surfaces  for  convenient  manipulation  and  geometry  modification. 
After  a  volume  grid  is  generated,  a  grid  pos^rocessing  program  POSTGRID  improves  the  grid 
quality  by  replacing  isolated  defects  in  the  volume  grid  using  local  remeshing.  The  USM3D  flow 
solver  was  developed  by  Neal  Frink  at  NASA  Langley  Research  Center.  The  methodology  is 
based  on  an  upwind,  tetrahedral  cell-centered  finite-volume  scheme  and  uses  a  Spalart-Allmaras 
one-equation  turbulence  model  coupled  with  a  wall  function  to  reduce  the  number  of  cells 
required  in  the  sublayer  region  of  the  boundary  layer.  Detailed  descriptions  and  results  from  the 
method  are  presented  in  Reference  (2). 

The  VGRID  grid  generation  system  was  distributed  by  NASA  Langley  Research  Center. 
The  distribution  uses  a  keying  system  whereby  each  software  executable  is  keyed  to  a  specific 
machine  and  ip  address.  The  VGRID  system  executables  are  only  available  for  Silicon  Graphics 
computer  systems.  The  USM3D  software  was  distributed  as  source  code  with  a  Makefile  valid 
for  several  computer  systems  (not  including  HP).  Very  minor  changes  were  required  to  the 
Makefile  to  compile  and  install  the  USM3D  software  on  the  MDC  Hewlett  Packard  (HP)  and 
government  supplied  CRAY  computer  systems. 

The  McDonnell  Aircraft  Computational  Grid  System  (MACGS)  is  a  versatile  interactive 
system  that  includes  capabilities  for  geometry  manipulation,  structured,  imstructured  and  chimera 
grid  generation,  and  boundary  condition  specification.  MACGS  consists  of  the  modules, 
ZONI3G,  GMAN,  and  GPRO.  ZONI3G  is  used  for  boundary  grid  generation,  GMAN  contains 
the  volume  grid  generation,  and  GPRO  performs  grid  manipulations  such  as  scaling 
transformations  etc.  The  NAvier  Stokes  Time  Dependent  (NASTD)  flow  solver  is  the  primary 
flow  solver  in  use  at  MDC/St.  Louis.  This  zonal  flow  solver  contains  capabilities  for  structured, 
unstructured  and  chimera  grids.  NASTD  and  MACGS  have  been  used  on  a  wide  range  of 
diverse  problems  at  MDC  and  at  several  government  sites.  These  tools  were  already  installed 
on  all  of  the  computer  systems  used  in  this  study. 

The  TETMESH  grid  generation  software  was  developed  by  Computational  Mechanics 
Company  under  contract  from  USAF  Wright  Laboratory.  TETMESH  uses  a  constrained 
Delaunay  procedure  for  generation  of  Euler  tetrahedral  grids  and  a  surface  extrusion  method  for 
the  viscous  layer  grid.  Cobalt  was  developed  at  the  USAF  Wright  Laboratory.  The  methodology 
is  based  on  an  upwind,  cell-centered  finite-volume  scheme  and  uses  a  Spalart-Allmaras  one- 
equation  turbulence  model.  Cobalt  can  use  grids  composed  of  cells  of  arbitrary  type  and  different 
cell  types  are  permitted  within  the  same  grid. 
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The  TETMESH  software  was  distributed  by  WL/FIMC  personnel  as  a  file  containing  the 
TETMESH  executable  and  several  utility  files.  The  software  only  runs  on  Silicon  Graphics 
computer  systems  and  was  installed  without  difficulty.  The  Cobalt  software  was  distributed  as 
sovirce  code  with  a  Makefile  appropriate  for  several  computer  systems  (not  including  HP). 
Several  minor  changes  were  required  to  the  Makefile  to  compile  and  install  the  software  on  the 
MDC  HP  and  government  supplied  CRAY  computer  systems. 
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3.  UNSTRUCTURED  GRID  ASSESSMENT:  FLAT  PLATE  FLOWFIELDS 


®xtensive  availability  of  test  and  analytical  data  for  flat  plate  flow  fields  makes  this 
evaluatmg  the  accuracy  of  flow  solver  algorithms.  In  addition,  the  simple 
e  of  the  flat  plate  geometry  reduces  grid  generation  and  solution  run  times.  These  properties 
^  ed  tK  to  generate  a  large  number  of  grids  and  solutions  in  a  relatively  short  amount  of  time. 
These  solutions  were  used  to  evaluate  the  relative  accuracy  and  to  identify  appropriate  grid 
spacmg  requirements  of  each  method  over  a  range  of  flow  field  conditions. 

viscous  layer  grid  requirements  of  each  flow  simulation  method  vary  due  to 
differences  m  numencal  algorithms.  Before  an  equitable  comparison  of  the  efficiency  and 
accuracy  of  each  method  can  be  made,  grid  spacing  requirements  must  be  established.  These 
requirements  were  identified  by  perfonning  a  grid  refinement  study  on  turbulent  and  laminar  flat 
p  ates  at  three  different  run  conditions.  Three  grid  spacing  parameters  were  determined  to  be 
cntacal  to  the  accuracy  of  the  flow  solution  for  each  method.  These  are  the  spacing  of  the  first 
^mt  normal  to  &e  wall,  the  number  of  points  across  the  boundary  layer,  and  the  surface  grid 
density.  The  grid  spacing  requirements  for  each  flow  simulation  method  were  determined  by 
sequentially  reducing  one  of  the  grid  spacing  parameters  and  generating  a  solution.  This  was 
continued  until  a  nearly  grid  independent  solution  was  obtained. 


Once  adequate  grid  spacings  requirements  were  identified,  a  relative  assessment  of  the 
accuracies  of  each  method  was  made.  Laminar  solutions  were  evaluated  to  study  the  diffusion 
properties  of  each  method  in  the  absence  of  turbulence  model  differences.  Differences  in  the 
turbulence  modeling  implementation  were  investigated  by  comparing  turbulent  flat  plate  solutions. 


3.1  Grid  Generation 

A  flat  plate  spanning  30  inches  long  and  3  inches  wide  was  selected  for  this  study.  The 
solution  domain  extends  5  inches  upstream  of  the  plate  leading  edge  and  protrudes  5  inches  into 
the  flowfield  normal  to  the  plate  surface  (See  Figure  2). 

The  flow  field  outside  the  boundary  layer  plays  a  minor  role  in  the  flow  solution  over 
a  flat  plate.  The  grid  topology  outside  the  viscous  layer,  therefore,  has  litfle  effect  on  the  flat 
plate  solution.  For  this  study,  grid  generation  was  simplified  by  extending  the  viscous  layer  grid 
topology  to  the  outer  boundary  so  that  a  single  grid  topology  exists  throughout  the  domain  {i.e. 
hexahedral,  tetrahedral  or  prismatic).  Flat  plate  hexahedral  grids  were  generated  naing  the 
MACGS  software.  Prismatic  grids  were  generated  by  extmding  a  surface  mesh  consisting  of 
triangular  faces  into  prisms.  Tetrahedral  flat  plate  grids  were  generated  by  subdividing  a 
hexahedral  grid  into  tetrahedra  (each  hexahedral  cell  was  subdivided  into  2  prismatic  cells,  which 
were  further  subdivided  into  3  tetrahedra  each  ).  A  total  of  31  flat  plate  grids  were  generated 
including  12  prismatic.  13  hexahedral,  and  6  tetrahedral  grids.  Three  measures  were  selected  to 
describe  the  grid  density.  These  are  spacing  of  the  first  point  normal  to  the  wall,  surface  grid 
density,  and  number  of  points  across  the  boundary  layer.  A  summary  of  the  grid  sizes,  and  grid 
density  for  each  grid  is  given  in  Table  1. 


5 


Prismatic 

2662 

696 

Prismatic 

696 

Prismatic 

696 

Prismatic 

696 

Prismatic 

696 

Prismatic 

696 

Prismatic 

696 

Grid 

Name 

Grid  Type 

SiorjEue 

^■:Faces  '  ■  • 

Noimal 

Pouts 

Nmiiberof 

Foutts 

Nomber  ; 

4^CeQs 

Ay  at  the 

waD 

n-6 

Hexahedral 

69x5 

31 

10,695 

8,160 

0.0001 

n-7 

Hexahedral 

69x5 

29 

10,005 

7,616 

0.00037 

n-8 

Hexahedral 

69x5 

27 

9315 

7,072 

0.00084 

n-9 

Hexahedral 

69x5 

24 

8,280 

6356 

0.0014 

n«10 

Hexahedral 

69x5 

59 

20355 

15,776 

0.00085 

n-11 

Hexahedral 

69x5 

59 

20355 

15,776 

0.0001 

n-12 

Hexahedral 

69x5 

41 

14,145 

10,880 

0.005 

n-13 

Hexahedral 

69x5 

41 

14,145 

10,880 

0.01 

n-14 

Hexahedral 

59x5 

41 

14,145 

10,880 

0.0028 

u-1 

Tetrahedral 

192 

15 

2205 

8064 

.0014 

u-2 

Tetrahedral 

192 

17 

2499 

9216 

.0004 

u-3 

Tetrahedral 

192 

13 

1911 

6912 

.01 

u-4 

Tetrahedral 

192 

12 

1764 

6336 

,03 

u-5 

Tetrahedral 

192 

12 

1764 

6336 

.02 

u-6 

Tetrahedral 

192 

13 

1911 

6912 

.0056 

3.2  Flow  Solutions 

Several  solutions  were  computed  on  the  grids  listed  in  Table  1  using  the  Cobalt,  NASTD, 
and  USM3D  flow  solvers.  These  methods  were  evaluated  across  a  range  of  flow  conditions. 
Laminar  flat  plate  solutions  were  generated  at  a  Mach  number  0.1  and  Reynolds  number  of  1.77 
million  at  the  plate  trailing  edge.  This  condition  provided  the  ability  to  evaluate  the  difliision 
properties  of  each  of  the  method  in  the  absence  of  turbulent  modeling  effects.  The  second  run 
condition  of  Mach  0.1  and  Reynolds  number  of  1.77  million  was  selected  to  provide  an 
assessment  of  a  turbulent  flowfield  at  a  low  to  moderate  Reynolds  number.  In  addition,  the 
relatively  low  Mach  number  of  0.1  provided  an  almost  incompressible  flow  field. 
Compressibility  effects  were  evaluated  by  computing  turbtilent  solutions  at  a  Mach  number  of 
0.7  and  Reynolds  number  of  12  million  at  the  plate  trailing  edge. 

A  total  of  57  flat  plate  solutions  were  computed  consisting  of  23  Cobalt,  24  NASTD,  and 
10  USM3D  solutions.  A  summary  of  the  run  conditions  for  each  solution  is  shown  in  Table  2. 
A  heavy  line  separates  solutions  from  each  run  condition  and  each  flow  solver.  In  addition  to 
providing  relative  assessments  of  accuracy,  each  solution  was  used  to  evaluate  one  or  more  grid 
meastires  as  part  of  the  grid  refinement  study.  The  grid  measures  that  were  assessed  with  a 
particular  solution  are  indicated  in  the  shaded  columns  with  a  +  for  spacing  of  the  first  grid  point 
normal  to  the  surface,  b  for  the  number  of  boundary  layer  points,  and  s  for  the  surface  grid 
density. 
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Table  2.  Flat  Plate  Run  Matrix 


NASTD 

NASTD 

NASTD 

NASTD 

NASTD 


1022 

8 

1.99 

13 

3.61 

9 

023 

22 

0.86 

15 

1.98 

13 

did  not  use 


+  b 

+  :  ■ 

+  s 
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An  attempt  was  made  to  minimize  the  differences  between  the  flow  solver  algorithms  by 
selecting  similar  algorithm  features  when  available.  For  instance,  all  solutions  were  generated 
using  a  second-order  accurate  (in  computational  space)  algorithm  and  the  Spalart/Allmaras 
turbulence  model  (for  turbtilent  computations).  Even  with  these  choices  several  differences 
remain  between  the  flow  solver  algorithms.  One  such  difference  is  in  the  flux  computations. 
Both  NASTD  and  USM3D  use  Roe’s  approximate  Reimann  solver  for  these  computations  while 
Cobalt  uses  the  approximate  Reimann  solution  method  due  to  Collela  in  combination  with  the 
iterative  method  of  Gottlieb  and  Groth.  There  are  also  differences  in  the  viscous  term 
computations  between  the  methods.  For  turbulent  flow  fields  USM3D  uses  a  wall  function  to 
model  the  near  wall  region.  The  time  integration  schemes  of  each  method  are  very  different. 
Cobalt  uses  an  explicit  Runge-Kutta  algorithm  while  NASTD  uses  an  approximate  factorization 
implicit  method  and  USM3D  uses  a  point  Jacobi  implicit  method.  The  spatial  discretization  is 
also  substantially  different  in  the  three  codes  with  Cobalt  and  USM3D  using  cell  centered 
methods  and  NASTD  using  a  vertex  based  scheme. 
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Solution  converge  of  all  solutions  was  determined  by  monitoring  the  integrated  drag  and 
skin  fnction.  ConvCTgence  data  for  each  solution  is  provided  in  Appendices  and  C. 

3.3  Laminar  Flat  Plate  Results 

Laminar  flat  plate  solutions  were  generated  at  a  Mach  nvunber  0.1  and  Reynolds  number 
of  1.77  million.  This  condition  provided  the  ability  to  evaluate  the  diffusion  properties  of  each 
of  the  method  in  the  absence  of  turbulence  modeling  effects.  Flow  solution  accuracy  was 
evaluated  by  comparing  the  predicted  velocity  profile  results  with  the  well  known  theoretical 
profile  of  Blasius.  The  skin  fnction  was  compared  with  the  theoretical  skin  friction  for  laminar 
flow  of  Cp  1.328  (ReJ’°-^ 

33.1  Cobalt  Results 

Cobalt  Laminar  flat  plate  solutions  were  run  on  grids  c-7,  c-8,  c-9,  c-11,  c-12,  and  c-13. 
The  solutions  on  grids  c-1 1  and  c-12  did  not  fiilly  converge  and  so  the  results  are  not  shown  in 
this  report.  All  of  these  grids  contain  the  same  number  of  surface  points.  Attempts  to  run  on 
grids  c-2  through  c-5  with  larger  surface  sizes  required  an  excessively  large  amount  of  time  to 
converge.  After  several  weeks  of  running  these  solutions  were  abandoned.  The  explicit  time 
integration  scheme  used  by  Cobalt  dictated  a  relatively  low  CFL  number  of  0.7.  As  a  result  the 
Cobalt  laminar  flat  plate  solutions  were  foimd  to  converge  very  slowly  requiring  upwards  of 
300,000  iterations. 

The  effect  of  wall  grid  spacing  on  solution  accuracy  is  illustrated  by  the  solutions  on  grids 
c-7,  c-9  and  c-13.  These  solutions  provide  wall  grid  spacings  of  y'^=0.39,  1.5  and  12 
respectively.  The  Cobalt  predicted  velocity  profile  across  the  boxmdary  layer  is  compared  with 
the  Blasius  profile  at  four  plate  locations  in  Figures  3  through  5.  The  profiles  from  grids  c-7, 
and  c-9  are  similar.  The  profiles  from  grids  c-9  and  c-13  differ  primarily  in  the  upstream  portion 
of  the  plate.  The  predicted  skin  fnction  from  these  solutions  is  compared  with  the  theoretical 
value  in  Figure  6.  The  wall  grid  spacing  was  found  to  have  little  effect  on  the  predicted  skin 
fnction  for  the  two  sizes  that  were  run.  The  solution  near  the  leading  edge  of  the  plate  is  slightly 
better  on  grid  c-9  (y*  of  1.5)  compared  to  grid  c-13  (y*  of  12).  Clearly  additional  data  is  needed 
to  determine  the  required  grid  spacing  for  these  run  conditions.  The  solution  on  grid  c-12  would 
have  provided  an  intermediate  data  point  for  the  wall  grid  spacing.  Unfortunately  this  solution 
did  not  fully  converge.  From  the  data  that  we  do  have  we  can  conclude  that  the  required  wall 
grid  spacing  for  laminar  solutions  is  somewhere  between  1.5  and  12. 
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Figure  3.  Velocity  profile,  Cobalt  laminar 
solution,  Mach  0.1,  Re=1.77xl0^  Grid  c-7. 


Figure  4.  Velocity  profile.  Cobalt  laminar 
solution,  Mach  0.1,  Re=1.77xl0«,  Grid  c-9. 


Figure  5.  Velocity  profile.  Cobalt  laminar 
solution,  Mach  0.1,  Re=1.77xl0®,  Grid  c-13. 


x{in) 


Figure  6.  Effect  of  wall  grid  spacing  on  skin  friction  coefficient. 

Cobalt  laminar  solutions,  Mach  0.1,  Re=1.77  million 

33.2  NASTD  Results 

NASTD  Laminar  flat  plate  solutions  were  run  on  grids  n-6,  n-7,  n-8,  n-9,  n-1 1,  n-12,  and 
n-13.  The  solutions  on  grids  n-7  and  n-8  were  virtually  identical  to  the  solution  on  grid  n-6  and 
so  the  results  are  not  shown.  Grids  n-9,  n-12  and  n-13  provided  a  comparison  of  solutions  with 
wall  grid  spacings  of  =  1.5,  6,  and  10  and  grids  n-6  and  n-11  provide  a  comparison  of  two 
different  distributions  of  points  across  the  boundary  layer.  All  of  these  grids  have  the  same 
surface  grid  size.  The  time  integration  was  performed  with  a  factored  implicit  scheme  permitting 
CFL  numbers  of  approximately  5.0.  The  NASTD  laminar  flat  plate  solutions  were  formd  to 
converge  slower  than  the  NASTD  turbulent  cases  requiring  between  30,000  and  80,000  iterations. 

The  effect  of  the  wall  grid  spacing  on  solution  accuracy  is  demonstrated  by  the  solutions 
on  grids  n-9,  n-12  and  n-13.  The  NASTD  predicted  velocity  profile  on  these  grids  at  four  plate 
locations  is  compared  with  the  Blasius  profile  in  Figures  8,  10,  and  11.  Predicted  skin  friction 
is  compared  with  the  theoretical  value  in  Figure  12.  Once  the  boundary  layer  was  developed, 
the  wall  grid  spacing  was  found  to  have  little  effect  on  the  predicted  skin  friction  for  the  three 
sizes  that  were  run.  The  solution  near  the  leading  edge  of  the  plate,  however,  was  predicted 
better  on  the  grid  with  a  y"^  of  1.5  compared  to  the  other  two  grids.  The  predicted  velocity 
profiles  also  agree  very  well  with  the  theoretical  Blasius  solution  once  the  boundary  layer  was 
developed. 

The  effect  of  the  number  of  points  across  the  boundary  layer  was  investigated  by 
comparing  the  solutions  on  grids  n-6  and  n-11.  These  grids  have  18  and  33  points  across  the 
viscous  layer.  The  predicted  velocity  profiles  using  these  grids  are  shown  in  Figures  7  and  9, 
and  Ihe  predicted  skin  friction  is  shown  in  Figure  13.  The  distribution  of  points  was  formd  to 
have  httle  effect  on  the  skin  friction  and  velocity  profile  predictions. 
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Figure  7.  Velocity  profile,  NASTD  laminar 
solution,  Mach  0.1,  Re=1.77xl0^  Grid  n-6. 
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Figure  8.  Velocity  profile,  NASTD  laminar 
solution,  Mach  0.1,  Re=1.77xl0^  Grid  n-9. 
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Figure  9.  Velocity  profile,  NASTD  laminar 
solution,  Mach  0.1,Re=L77xl0®,  Grid  n-11. 
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Figure  10.  Velocity  profile,  NASTD  laminar 
solution,  Mach  0.1,  Re=l. 77x10®,  Grid  n- 
12. 
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Fi^e  11.  Velocity  profile,  NASTD 
laminar  solution,  Mach  0,1, 
Re=1.77xlO®,Grid  n-13 
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12.  Effect  of  wall  grid  spacing  on  skin  friction  coefficient,  NASTD  laminar  solutions, 

Mach  0.1,  Re=1.77  million. 


X  (in) 

Figure  13.  Effect  of  number  of  points  in  botmdary  layer  on  skin  friction  coefficient,  NASTD 

laminar  solutions,  Mach  0.1,  Re=1.77  milimn 
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333  USM3D  Results 

USM3D  laminar  flat  plate  solutions  were  run  on  grids  u-3  and  u-6  providing  a  comparison 
of  solutions  with  wall  grid  spacings  of  y"  =  6.25  and  12.22.  Both  solutions  have  the  same 
surface  grid  size.  The  time  integration  was  performed  with  a  point  Jacobi  implicit  scheme 
permittmg  CFL  numbers  of  approximately  200.  The  solutions  on  grid  u-3  and  u-6  required 
20,000  and  12,000  iterations  to  converge. 

USM3D  predicted  velocity  profiles  using  these  grids  are  shown  in  Figures  14  and  15. 
The  skin  friction  predictions  are  compared  with  the  theoretical  value  in  Figure  16.  Once  the 
boundary  layer  is  developed,  the  wall  grid  spacing  has  little  effect  on  the  predicted  skin  friction 
for  the  two  wall  spacing  sizes  that  were  run.  The  solution  near  the  leading  edge  of  the  plate, 
however,  is  predicted  better  on  the  grid  with  tighter  wall  spacing.  The  predicted  velocity  profiles 
agree  very  well  with  theoretical  Blasius  solution  once  the  boundary  layer  is  developed.  The 
profiles  near  the  leading  edge  do  not  agree  well  with  the  Blasius  solution,  probably  to  poor 
grid  resolution. 


y  sqrt(uinf^x) 

Figure  14.  Velocity  profile,  USM3D  lamiTiar 
solution,  Mach  0.1,  Re=  1.77x10®,  Grid  u-3. 
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Figure  15.  Velocity  profile,  USM3D  laminar 
solution,  Mach  0.1,  Re=1.77xl0®,  Grid  u-6. 
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Figure  16.  Effect  of  wall  grid  spacing  on  skin  friction  coefficient, 
USM3D  laminar  solutions,  Mach  0.1,  Re=1.77  million. 


3.4  Turbulent  Flat  Plate  Results 

Tiirbulent  flat  plate  solutions  were  computed  at  two  run  conditions  with  all  three  of  the 
flow  solution  methods.  Nearly  incompressible  solutions  were  computed  at  a  run  condition  of 
Mach  0. 1 ,  Reynolds  number  of  1 .77  million  while  the  effects  of  compressibility  were  investigated 
at  run  conditions  of  Mach  0.7  and  Reynolds  number  of  12  million.  All  of  the  solutions  were 
computed  using  the  Spalart/Alhnaras  turbulence  model.  Flow  solution  accuracy  was  evaltiated 
by  comparing  &e  predicted  velocity  profile  results  with  the  theoretical  law  of  the  wall  velocity 
profile.  The  skin  ffiction  was  compared  with  the  theoretical  skin  friction  for  fully  turbulent  flow 
of  Cf=0.0583  (ReJ*’'^  which  is  based  on  Prandtl’s  power  law  assximption.  For  the  hi^  Reynolds 
number  case,  the  skin  friction  is  also  compared  against  the  relationship  Cp0.455/(ln^(0.06  ReJ) 
which  is  taken  from  Reference  (3). 

3.4.1  Cobalt  Results 

After  running  the  Cobalt  flat  plate  solutions  to  convergence,  a  change  in  the  Cobalt 
implementation  of  the  Spalart/Alhnaras  turbulence  model  was  made  in  mid  October  1996.  Rather 
than  rerun  all  of  the  turbulent  Cobalt  solutions  from  scratch,  these  solutions  were  restarted  from 
their  previous  states.  The  number  of  iterations  shown  in  Table  2  includes  only  the  iterations 
required  to  reach  convergence  before  the  software  change.  The  software  change  resulted  in  a 
slight  change  in  the  converged  value  of  integrated  drag  as  shown  in  the  drag  convergence  plots 
in  Appendix  A. 

Turbulent  flat  plate  solutions  at  the  first  run  condition  of  Mach  0.1  were  generated  on 
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grids  u-1,  u-2,  c-6,  c-7,  c-8,  c-9,  c-11,  and  c-14.  Solutions  on  grids  c-3  and  c-4  were  also  run 
but  were  not  restarted  after  the  turbulaice  model  modification  was  made  to  Cobalt.  The  solution 
on  grid  c-11  would  not  converge  and  was  discontinued  after  163,000  iterations.  The  reason  for 
this  lack  of  convergence  is  not  known  although  it  may  be  grid  related  since  the  laminar  solution 
on  this  grid  also  failed  to  converge.  These  solutions  converged  relatively  slowly  requiring 
approximately  150,000  iterations  to  reach  a  converged  solution. 

The  solutions  on  grids  c-7,  c-8,  c-9,  and  c-14  have  wall  grid  spacings  of  y*=0.58,  y'^=1.32, 
y'^=2.24,  and  y^=4.48.  Predicted  velocity  profiles  across  the  boimdary  layer  at  four  plate  locations 
are  compared  with  the  law  of  the  wall  in  Figures  1 7  through  2 1 .  The  predicted  skin  fiiction  from 
these  solutions  is  compared  with  the  theoretical  value  in  Figure  22.  The  skin  fnction  changes 
as  the  wall  grid  spacing  is  refined  until  a  spacing  of  y'^=1.32  is  reached.  The  predicted  skin 
fnction  and  velocity  profiles  at  wall  spacings  of  y*=0.58  and  1.32  are  nearly  identical.  We  were 
unable  to  determine  the  effect  of  grid  density  across  the  boundary  layer  due  to  lack  of 
convergence  of  the  solution  on  grid  c-11. 

The  effect  of  grid  cell  shape  on  solution  accuracy  was  evaluated  by  comparing  Cobalt 
solutions  on  prismatic  cell  grids  c-7  and  c-9  with  solutions  generated  on  tetrahedral  cell  grids  u-2 
and  u-1.  The  wall  grid  spacing  for  grids  c-7  and  u-2  are  0.00037  and  0.0004  which  results  in 
a  y"^  of  approximately  0.6  at  these  run  conditions.  The  wall  grid  spacing  for  grids  c-9  and  u-2 
is  0.0014  resulting  in  a  y*  of  approximately  2.24.  The  predicted  skin  fiiction  for  these  solutions 
is  compared  with  the  theoretical  value  in  Figure  23.  The  solutions  on  grids  c-9  and  u-1  are 
nearly  identical.  However,  the  solution  on  tetrahedral  grid  u-2  is  closer  to  the  theoretical  value 
than  the  solution  on  prismatic  grid  c-7.  This  indicates  that  cell  shape  does  influence  solution 
accuracy  in  some  instances. 


Figure  17.  Velocity  profile.  Cobalt  solution.  Figure  18.  Velocity  profile.  Cobalt  solution, 

Mach  0.1,  Re=1.77  million.  Grid  c-6.  Mach  0.1,  Re=1.77  million.  Grid  c-7. 
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Figure  19.  Velocity  profile,  Cobalt  solution, 
Mach  0.1,  Re=1.77  million,  Grrid  c-8. 


Figure  21.  Velocity  profile.  Cobalt  solution, 
Mach  0.1,  Re=1.77  million.  Grid  c-14. 


Mach  0.1,  Re=1.77  million.  Grid  c-9. 
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solutions,  Mach  0.1,  Re=1.77  million 


Figure  23.  Effect  of  cell  shape  on  predicted  skin  friction  coefficient,Cobalt  turbulent 

solutions,  Mach  0.1,  Re=1.77  million 
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Cobalt  turbulent  flat  plate  solutions  were  run  at  a  second  run  condition  of  Mach  0.7,  and 
Rejmolds  number  of  12  million.  At  these  run  conditions,  solutions  were  generated  on  grids  u-2, 
c-6,  c-7,  c-8,  c-9,  and  c-11.  Solutions  on  grids  c-3,  c-4,  and  c-5  were  also  run  but  were  not 
restarted  after  the  turbulence  model  modification  was  made  to  Cobalt.  These  solutions  converged 
much  faster  than  the  low  Reynolds  number  solutions  requiring  approximately  60,000  iterations. 

The  solutions  on  grids  c-6,  c-7,  c-8,  and  c-9  have  wall  grid  spacings  of  y*=0.8,  y‘^=3.0, 
y*=7.7,  and  y'*’=12.8.  These  solutions  were  used  to  evaliiate  the  effect  of  wall  grid  spacing  on 
solution  accuracy.  Velocity  profiles  across  the  viscous  layer  from  the  Cobalt  solutions  are  shown 
in  Figures  24-27.  The  predicted  skin  fnction  from  these  solutions  is  compared  with  the 
theoretical  value  in  Figure  29.  The  solutions  do  not  appear  to  change  at  wall  grid  spacings  below 
y*=3.0. 

The  solutions  on  grids  c-6  and  c-11  were  used  to  evaluate  the  effect  of  viscous  layer  grid 
resolution  on  solution  accuracy.  These  grids  contain  17  and  40  points  across  the  viscous  layer 
{y*  <  800  )  respectively.  Predicted  velocity  profile  and  skin  fiiction  from  the  two  solutions  is 
shown  in  Figures  24,28  and  30.  Comparison  of  the  skin  friction  and  velocity  profiles  from  the 
two  solutions  indicates  that  adequate  viscous  layer  grid  resolution  is  achieved  with  17  grid  points. 

The  effect  of  cell  shape  on  solution  accuracy  was  evaluated  by  comparing  the  solutions 
on  grids  c-7  and  u-2.  These  grids  have  a  wall  spacing  of  .0004  which  results  in  a  y"^  of 
approximately  3.0  at  these  run  conditions.  The  predicted  skin  friction  for  these  two  solutions  is 
compared  with  the  theoretical  value  in  Figure  31.  Once  again  the  cell  shape  does  not  have  a 
large  influence  on  the  solution  accuracy. 
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Figure  24.  Velocity  profile,  Cobalt  solution,  Figure  27.  Velocity  profile.  Cobalt  solution, 
Mach  0.7,  Re=12  million.  Grid  c-6.  Mach  0.7,  Re=12  million.  Grid  c-9. 


Figure  25.  Velocity  profile.  Cobalt  solution,  Figure  28.  Velocity  profile,  Cobalt  solution, 
Mach  0.7,  Re=12  million.  Grid  c-7.  Mach  0.7,  Re=12  million.  Grid  c-11. 
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Figure  26.  Velocity  profile.  Cobalt  solution, 
Mach  0.7,  Re=12  milhon.  Grid  c-8. 
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29.  Effect  of  wall  grid  spacing  on  predicted  skin  friction  coefficient,  Cobalt  turbulent 
solutions,  Mach  0.7,  Re=12  million. 


Figure  30.  Effect  of  number  of  boundary  layer  points  on  predicted  skin  fiiction  coefficient. 
Cobalt  turbulent  solutions,  Mach  0.7,  Re=12  million. 
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Figure  3 1 .  Effect  of  cell  shape  on  predicted  skin  friction  coefficient.  Cobalt  turbulent 

solutions,  Mach  0.7,  Re=12  million. 


3.4.2  NASTD  Results 

NASTD  turbulent  flat  plate  solutions  were  generated  at  the  first  run  condition  of  Mach 
0.1,  Re5molds  number  of  1.77  million  on  grids  n-4,  n-6,  n-7,  n-8,  n-9,  n-11,  and  n-14.  These 
solutions  were  run  at  a  CFL  number  of  5.0  and  required  between  10,000  and  40,000  iterations 
to  reach  convergence. 

The  effect  of  wall  grid  spacing  on  solution  accuracy  was  evaluated  by  computing  solutions 
on  grids  n-7,  n-8,  n-9,  and  n-14.  At  this  run  condition  these  grids  provide  wall  grid  spacings  of 
y'^=0.9,  y^=2.0,  y'^=4.0,  and  y‘^=8.0.  The  NASTD  predicted  velocity  profile  across  the  viscous 
layer  for  these  grids  is  compared  with  die  law  of  the  wall  in  Figures  33  through  35  and  37.  The 
predicted  skin  fiiction  from  these  solutions  is  compared  with  ^e  theoretical  value  in  Figure  38. 
The  velocity  profiles  and  skin  friction  reach  a  nearly  grid  independent  solution  at  a  wall  spacing 
of  y'^=2.0. 

The  effect  of  viscous  layer  grid  resolution  was  evaluated  using  grids  n-6  and  n-11. 
Solutions  on  these  grids  contain  22  and  40  points  across  the  viscous  layer  (y"^  <  800)  respectively. 
The  predicted  velocity  profile  across  the  viscous  layer  for  these  grids  is  shown  in  Figures  32  and 
36.  The  predicted  skin  friction  is  compared  with  the  theoretical  value  in  Figure  39.  The  two 
solutions  are  very  similar  indicating  that  adequate  grid  resolution  is  achieved  with  22  points. 

The  effect  of  surface  grid  resolution  was  evaluated  by  comparing  NASTD  solutions  on 
grids  n-4  and  n-8.  These  grids  have  surface  grid  sizes  of  118x5  and  69x5  respectively.  The  skin 
friction  predicted  using  grids  n-4  and  n-8  is  compared  in  Figure  40.  The  similarity  of  the  two 
solutions  indicates  that  the  coarser  surface  grid  provides  adequate  resolution. 
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Figure  32.  Velocity  profile,  NASTD 
solution,  Mach  0.1,  Re=1.77xl0^  Grid  n-6. 


Figure  3 5. Velocity  profile,  NASTD  solution, 
Mach  0.1,  Re=1.77xl0^  Grid  n-9. 


Figure  33.  Velocity  profile,NASTD  solution, 
Mach  0.1,  Re=l. 77x10®,  Grid  n-7. 


Figure  34.  Velocity  profile,  NASTD  solution, 
Mach  0.1,  Re=l. 77x10®,  Grid  n-8. 


Figure  36.  Velocity  profile,  NASTD 
solution,  Mach  0.1,  Re=l. 77x10®,  Grid  n-1 1. 


Figure  37.  Velocity  profile,  NASTD 
solution,  Mach  0.1,  Re=l  .77x10®,  Grid  n-14. 
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e  38.  Effect  of  wall  grid  spacing  on  skin  Mction  coefficient,  NASTD  turbtilent  solutions 

Mach  0.1,  Re=1.77  million. 
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Figure  39.  Effect  of  number  of  boundary  layer  points  on  skin  friction  coefficient,  NASTD 
turbulent  solutions,  Mach  0.1,  Re=1.77  million. 


solutions,  Mach  0.1,  Re=1.77  million. 

NASTD  turbxilent  flat  plate  solutions  were  run  at  the  second  run  condition  of  Mach  0.7, 
and  Reynolds  number  of  12  million.  At  these  run  conditions,  solutions  were  generated  on  grids 
n-2,  n-6,  n-7,  n-8,  n-9,  and  n-11.  These  solutions  converged  relatively  quickly  requiring 
approximately  4,000  iterations. 

The  required  wall  grid  spacing  was  evaluated  by  computing  solutions  on  grids  n-6,  n-7, 
n-8,  and  n-9.  These  grids  have  wall  spacings  of  y'"=1.3,  y*=4.9,  y'^=12.4,  and  y*=23.6.  The 
predicted  velocity  profile  for  these  solutions  is  shown  in  Figures  41  through  44.  The  skin  friction 
is  compared  with  the  theoretical  value  in  Figure  46.  The  predicted  skin  friction  is  different  on 
all  four  grids.  This  indicates  that  a  y"*^  value  below  4.9  is  necessary  to  reach  a  nearly  grid 
independent  solution. 

The  effect  of  viscous  layer  grid  resolution  was  evaluated  with  the  solutions  on  grids  n-6 
and  n-11.  These  grids  contain  17  and  30  points  across  the  viscous  layer  (y"^  <  800 )  respectively 
at  the  given  run  condition.  Velocity  profile  predictions  from  NASTD  solutions  on  these  grids 
are  shown  in  Figures  41  and  45  and  the  skin  friction  is  shown  in  Figure  47.  The  skin  friction 
and  velocity  profiles  from  the  two  solutions  are  very  similar  indicating  that  adequate  grid 
resolution  is  achieved  with  17  points.  The  effect  of  surface  mesh  resolution  was  evaluated  by 
comparing  solutions  on  grids  n-2  and  n-6.  These  grids  contain  surface  grid  sizes  of  1 1 8x5  and 
69x5  respectively.  Skin  friction  predictions  for  these  two  solutions  are  shown  in  Figure  48.  The 
siimlarity  of  these  two  solutions  indicates  the  coarser  surface  grid  provides  adequate  grid 
resolution. 
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Figure  41.  Velocity  profile,  NASTD 
solution,  Mach  0.7,  Re=12  million,  Gridn-6. 


Figure  42.  Velocity  profile,  NASTD 
solution,Mach  0.7,  Re=12  million.  Grid  n-7 


solution,  Mach  0.7,  Re=12  million.  Grid  n-8. 


Figure  44.  Velocity  profile,  NASTD 
solution,  Mach  0.7,  Re=12  million.  Grid  n-9 


Mach  0.7,  Re=12  million.  Grid  n-11. 
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Figure  46.  Effect  of  wall  grid  spacing  on  skin  friction  coefficient, NASTD  turbulent  solutions, 

Mach  0.7,  Re=12  million. 


Figure  47.  Effect  of  number  of  boundary  layer  points  on  skin  friction  coefficient,  NASTD 

turbulent  solutions,  Mach  0.7,  Re=12  million. 
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Figure  48.  Effect  of  surface  grid  resolution  on  skin  friction  coefficient,  NASTD  turbulent 

solutions,  Mach  0.7,  Re=12  million. 


3,4,3  USM3D  Results 

USM3D  turbulent  flat  plate  solutions  were  generated  at  the  first  run  condition  of  Mach 
0.1,  Reynolds  number  of  1.77  million  on  grids  u-2,  u-3,  u-5,  and  u-4.  These  solutions  required 
between  5,000  and  1 0,000  iterations  to  reach  solution  convergence. 

The  required  grid  spacing  at  the  wall  was  evaluated  by  comparing  solutions  on  grids  u-2, 
u-3,  u-5,  and  u-4.  These  grids  have  wall  grid  spacings  of  y‘'=l.l,  y'"=25.4,  y'^=49.0,  and  y^=73.8 
respectively.  Predicted  velocity  profiles  across  the  viscous  layer  are  compared  with  the  law  of 
the  wall  at  four  locations  on  the  plate  in  Figures  49  through  52.  The  predicted  skin  friction  from 
each  solution  is  compared  with  the  theoretical  value  in  Figure  53.  For  y"*"  values  above  49.0,  the 
predicted  skin  friction  compares  very  well  with  the  theoretical  values.  At  these  large  wall  grid 
spacings,  the  skin  friction  is  primarily  governed  by  the  wall  function.  The  predicted  skin  friction 
on  grids  u-2  and  u-3,  however,  does  not  agree  as  well.  At  these  tighter  wall  grid  spacings  the 
wall  function  has  a  smaller  influence  on  the  near  wall  solution.  The  poor  solution  at  the  tighter 
grid  resolution  seems  to  indicate  lack  of  convergence  of  the  solution  or  improper  treatment  of 
the  solution  computation  on  highly  stretched  grid  cells.  Reasons  for  this  discrepancy  were  not 
investigated  beyond  running  the  solution  at  a  lower  CFL  number  which  had  no  effect. 

USM3D  turbulent  flat  plate  solutions  were  generated  at  the  second  run  condition  of  Mach 
0.7,  and  Reynolds  ntunber  of  12  million  on  grids  u-1,  u-6,  and  u-3.  These  solutions  converged 
relatively  quickly  requiring  approximately  4,000  iterations  at  a  CFL  number  of  200. 

Wall  grid  spacing  sensitivity  was  evaluated  by  comparing  the  three  solutions.  These 
solutions  have  wall  grid  spacings  of  y'^=18.2,  y'^=68.9,  and  y'^=122.4  respectively.  The  predicted 
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velocity  profile  is  shown  in  Figures  54  through  56  and  the  skin  friction  is  shown  in  Figure  57. 
Once  again  the  skin  friction  is  in  excellent  agreement  with  the  theoretical  values  for  grid  spacings 
above  y'^=49.0.  However,  at  the  smallest  wall  grid  spacing  where  the  wall  frmction  is  losing  its 
mfluence,  the  skin  friction  agreement  is  worse. 


Figure  49  .Velocity  profile,  USM3D  turbulent 
solution,  Mach  0.1,  Re=l. 77x10®,  Grid  u-2. 


Figure  5 1  .Velocity  profile,  USM3D  turbulent 
solution,  Mach  0.1,  Re=l. 77x10®,  Grid  u-5. 


Figure  50.  Velocity  profile,USM3D  turbulent 
solution,  Mach  0.1,  Re=l. 77x10®,  Grid  u-3. 


Figure  52.  Velocity  profile,USM3D  turbulent 
solution,Mach  0.1,  Re=1.77xl0®,  Grid  u-4. 
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Figure  53.  Effect  of  wall  grid  spacing  on  skin  friction  coefficient,  USM3D  turbulent  solutions, 

Mach  0.1,  Re=L77  million. 


Figure  54.Velocity  profile,  USM3D  turbulent  Figure  55.Velocity  profile,  USM3D  turbiilent 

solution,  Mach  0.7,  Re=12xl0^  Grid  u-1.  solution,  Mach  0.7,  Re=12xl0^  Grid  u-6. 


31 


Figure  56.  Velocity  profile,  USM3D  turbulent 
solution,  Mach  0.7,  Re=12xl0^  Grid  u-3. 


Figure  57.  Effect  of  wall  grid  spacing  on  skin  fnction  coefficient,USM3D  turbulent  solutions, 

Mach  0.7,  Re=12  million. 


3.5  Unstructured  Grid  Requirements 

l^e  grid  requirements  identified  in  the  previous  section  are  summarized  below.  Most  of 
the  requirements  that  are  given  in  this  section  are  given  as  ranges  rather  than  specific  numbers. 
These  requirements  are  meant  to  be  used  as  guidelines  in  the  development  of  future  grids.  Since 
these  grids  will  most  likely  not  be  over  flat  plates,  the  given  requirements  should  only  be  used 
as  a  reference  point.  These  requirements  are  also  only  valid  for  the  specific  grid  cell  shape  that 
was  tested  with  each  flow  solver  (i.e.  prism  for  Cobalt,  hexahedral  for  NASTD,  tetrahedral  for 
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USM3D). 

A  siunmary  of  the  grid  reqxxirements  identified  for  each  flow  solver  is  shown  in  Table  3. 
Two  grid  parameters  are  given  for  each  case.  These  are  the  distance  between  the  wall  and  the 
first  point  normal  to  the  surface  (wall  spacing),  and  the  number  of  points  across  the  boundary 
layer  identified  as  the  number  of  points  below  y^=800.  A  blank  entry  in  Table  3  indicates  the 
parameter  was  not  determined  for  that  case. 


Table  3.  Flat  Plate  Grid  Requirements 


■ 

1.77 

laminar 

12.2  <y"< 

1.77 

>22 

12.0 

turbulent 

3.0  <y"<  7.7 

>  17 

1.77 

laminar 

1.77 

turbulent 

2.0  <  y-K  4.0 

>22 

>272 

turbulent 

>17 

>272 

1.77 

laminar 

12.2  <  y^< 

>-5 - 

1.77 

>  5 

12.0 

122.4  <  y"< 

>5 

3.6  Flow  Solver  Comparisons 

A  comparison  of  the  accuracy  and  efficiency  of  each  of  the  three  imstructured  grid 
methods  is  presented  in  this  section.  The  comparison  is  based  on  the  data  generated  in  the  grid 
refinement  study  presented  above. 

3.6.1  Accuracy 

The  different  grid  topologies  used  with  each  flow  solver  precludes  an  equitable  one-to-one 
comparison  of  flow  solver  accuracy  independent  of  the  grid.  However,  an  indication  of  the 
relative  accuracy  of  each  flow  solver/  grid  combination  can  be  made. 

The  results  obtained  at  the  laminar  run  condition  of  Mach  0.1,  Re=1.77  provides  the 
opportunity  to  compare  flow  solvers  in  the  absence  of  turbulence  model  effects.  The  predicted 
skin  friction  from  each  of  the  methods  is  compared  with  the  Blasius  solution  in  Figure  58.  The 
Cobalt  and  NASTD  results  are  both  very  close  to  the  Blasius  solution  with  the  USM3D  solution 
being  the  worst  of  the  three.  The  Cobalt  solution  does  a  slightly  bettCT  job  of  predicting  the  skin 
friction  near  the  leading  edge  of  the  plate.  The  predicted  velocity  profiles  near  the  trailing  edge 
of  the  plate  are  shown  in  Figure  59.  All  three  methods  do  a  good  job  of  matching  the  Blasius 
profile. 
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Figure  59.  Comparison  of  Predicted  Laminar  Velocity  Profile,  Mach  0.1,  Re— 1.77  Million. 
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The  second  set  of  run  conditions,  Mach  0.1,  Re=1.77  million,  provides  a  comparison  of 
the  three  methods  in  a  turbulent  nearly  incompressible  flow.  In  Figure  60  the  predicted  skin 
friction  is  compared  with  the  theoretical  skin  friction  for  fully  turbulent  flow.  The  USM3D 
solution  with  the  waU  function  does  a  very  good  job  of  matching  the  theoretical  value.  The 
Cobalt  and  NASTD  results  differ  from  the  theoretical  value  near  the  leading  edge  and  approach 
the  theoretical  value  towards  the  trailing  edge.  Velocity  profiles  from  the  three  methods  are 
compared  with  the  law  of  the  wall  profile  in  Figure  61.  All  three  methods  match  the  theoretical 
profile  reasonably  well. 

Cobalt  and  USM3D  solutions  were  generated  on  grid  u-2  at  Mach  0.1,  Re=1.77  million. 
These  solutions  provide  an  equitable  comparison  of  flow  solver  accuracy  since  any  grid 
dependencies  have  been  removed.  The  predicted  skin  friction  from  the  two  solutions  is  shown 
in  Figure  62.  Neither  method  matches  &e  theoretical  values.  However,  the  Cobalt  results  are 
smooth  and  have  the  correct  trends  whereas  the  USM3D  results  are  very  erratic. 


Figure  60.  Comparison  of  Predicted  Turbulent  Flat  Plate  Skin  Friction,  Mach  0.1,  Re=1.77xl0® 
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Figure  61.  Comparison  of  Predicted  Tiirbiilent  Velocity  Profile,  Mach  0.1,  Re=1.77x  10* 
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The  final  run  condition  of  Mach  0.7,  Re— 12  million,  provides  a  comparison  of  the  three 
methods  in  a  high  Reynolds  number  turbulent  flow.  The  skin  fiiction  from  the  three  methods 
IS  compared  with  the  theoretical  value  in  Figure  63.  Once  again  the  USM3D  results  with  the 
wall  fimction  compare  very  well  with  the  theoretical  value.  The  NASTD  and  Cobalt  results  also 
compare  reasonably  well.  The  predicted  velocity  profiles  are  compared  with  the  law  of  the  wall 
profile  in  Figure  64.  All  three  methods  provide  similar  results. 


Figure  63.  Comparison  of  Predicted  Turbulent  Skin  Friction,  Mach  0.7,  Re=  12x10^ 


Figure  64.  Comparison  of  Predicted  Turbulent  Velocity  Profile,  Mach  0.7,  Re=12xl0^ 
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3.6.2  Efficiency 

It  is  difficult  to  compare  the  efficiencies  of  the  flow  solvers  based  on  the  solutions 
computed  in  this  study  due  to  differences  in  grid  sizes.  However,  an  estimate  of  the  comparable 
run  times  for  the  different  solvers  can  be  made  by  comparing  the  time  per  iteration  per  cell  and 
the  total  number  of  iterations  required  to  reach  convergence.  This  comparison  does  not  include 
penalties  associated  with  the  gnd  cell  topology.  For  instance,  a  tetrahedral  cell  grid  will  contain 
a  larger  number  of  cells  in  the  boundary  layer  than  a  hexahedral  cell  grid.  A  discussion  of  these 
penalties  is  given  below. 

The  run  times  per  iteration  per  cell  for  each  of  the  flow  solvers  is  given  in  Table  2. 
These  numbers  were  obtained  by  averaging  run  times  for  several  solutions  on  an  HP- 735 
workstation.  If  we  multiply  these  numbers  by  the  number  of  iterations  required  to  reach 
convergence  we  obtain  an  overall  solution  time  per  cell  for  each  flow  solver.  Comparing  the 
numbers  for  the  turbulent  flat  plate  solutions  at  a  Reynolds  number  of  1 .77  million  we  find  that 
the  Cobalt,  NASTD,  and  USM3D  solvers  take  roughly  47,  24,  and  20  seconds  per  solution  per 
cell.  These  times  were  obtained  using  an  average  number  of  iterations  per  solution  of  160,000 
(Cobalt),  25,000  (NASTD),  and  12,000  (USM3D)  iterations.  For  an  equal  size  grid,  there  is  a 
ratio  of  2:1  between  Cobalt  and  NASTD  or  USM3D  solution  times.  This  difference  is  most 
likely  due  to  the  explicit  time  integration  scheme  of  Cobalt  versus  the  implicit  schemes  in 
NASTD  and  USM3D.  When  taking  advantage  of  the  wall  function,  (grids  with  y"^  25),  USM3D 
was  able  to  reduce  the  run  time  by  another  factor  of  2. 

There  was  a  large  variation  in  the  memory  requirements  for  each  flow  solver.  Once  again 
the  comparisons  will  be  made  on  a  per  cell  basis  ignoring  any  penalties  associated  with  the  grid 
cell  topology.  For  the  laminar  flat  plate  solutions  the  Cobalt,  NASTD  and  USM3D  flow  solvers 
each  required  180,  95  and  253  single  precision  words  of  memory  per  cell.  For  the  turbulent  flat 
plate  solutions  these  solvers  required  185,  97  and  256  single  precision  words  of  memory  per  cell. 
The  NASTD  memory  requirements  are  much  lower  than  for  Cobalt  and  USM3D.  This  is 
primarily  due  to  the  lower  memory  requirements  of  a  structured  grid  data  set.  USM3D  has  the 
largest  memory  requirements  most  likely  due  to  the  implicit  time  integration  scheme. 

The  grid  cell  topology  can  play  a  large  role  in  the  efficiency  of  the  flow  solution.  This 
is  primarily  dependent  on  the  number  of  cells  required  to  resolve  the  viscous  layer  and  the 
increase  in  work  per  cell  to  compute  the  fluxes  across  each  cell  face.  Intuitively  it  would  seem 
that  a  tetrahedral  grid  will  require  the  largest  number  of  cells  to  achieve  a  given  level  of 
resolution  while  prismatic  and  hexahedral  cells  should  allow  the  most  efficient  modeling  of  the 
viscous  layer.  In  this  study,  however,  grids  containing  prismatic  cells  were  much  larger  than  the 
corresponding  tetrahedral  or  hexahedral  grids.  This  was  driven  by  the  attempt  to  maintain 
isotropic  triangular  faces  near  the  leading  edge  of  the  plate  where  the  streamwise  spacing  was 
reduced.  The  hexahedral  and  anisotropic  tetrahedral  grids  were  able  to  resolve  this  streamwise 
spacing  without  increasing  the  resolution  in  the  crosswise  direction.  Ftirther  study  of  the  effect 
of  grid  cell  type  on  grid  size  is  needed. 
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4.  UNSTRUCTURED  GRID  ASSESSMENT:  ONERA  WING 


The  ONERA  M6  wing  has  been  used  extensively  throughout  the  CFD  community  for  the 
evaluation  and  validation  of  CFD  codes.  The  availability  of  an  excellent  set  of  test  data  for  this 
geometry  makes  it  well  suited  for  evaluating  CFD  codes  across  a  range  of  flow  conditions  (4). 
Viscous  unstructured  grid  solutions  were  computed  at  two  run  conditions  with  each  of  the  three 
unstractured  grid  methods.  The  flowfield  at  the  first  run  condition  of  Mach  0.84,  a=0.04, 
Re=l  1.7  million  is  subsonic  and  remains  attached  across  the  entire  wing  span.  The  flowfield  at 
the  second  condition  of  Mach  0.84,  a=5.06,  Re=11.7  million  is  transonic  with  a  shock  induced 
separation  over  the  outer  one  third  of  the  wing.  A  total  of  1 1  solutions  were  computed  on  7 
grids. 


4.1  Grid  Generation 

Grids  were  generated  about  the  ONERA  wing  using  each  of  the  three  grid  generator 
packages.  Mixed  cell  grids  containing  prismatic  and  tetrahedral  cells  were  generated  within 
TETMESH.  Zonal  hybrid  grids  containing  hexahedral  and  tetrahedral  grids  were  generated  in 
MACGS.  Anisotropic  tetrahedral  grids  were  generated  in  VGRED.  A  summary  of  the  seven 
grids  generated  for  this  study  is  given  in  Table  4.  Each  grid  name  provides  three  pieces  of 
information,  the  first  letter  indicates  the  grid  generator  that  was  used  to  create  the  grid,  the 
second  letter  indicates  whether  the  grid  is  appropriate  for  euler  (inviscid)  or  viscous  analysis,  and 
the  last  number  is  used  to  distinguish  between  multiple  grids  containing  the  same  first  two  letters. 
For  instance  grid  tv2  is  the  second  viscous  packed  grid  generated  in  TETMESH. 


Table  4.  ONERA  Wing  Grid  Summary 


Grid 

Grid  Generator 

Sur&ce  Size 

Number  of  Cells 

Number  of 
Pmnts 

nel 

15,279 

231,507 

42,410 

N/A 

tel 

TETMESH 

11,568 

163,930 

30,308 

N/A 

te2 

TETMESH 

32,510 

220,251 

43,807 

N/A 

tvl  ! 

TETMESH 

20,160" 

422,614" 

175,522" 

0.000015 

mvl 

3,264 

366,273 

138,148 

0.000031 

v\'l 

VGRID 

6,483 

356,472 

62,883 

0.000203 

v\'2 

VGRID 

6,483 

414,038 

96,116 

0.000060 

-  sizes  are  tor  tuU  wmg  gnd  (no  synunetry  plane) 


4.2  Cobalt/TETMESH  Results 

A  total  of  5  Cobalt  solutions  were  computed  on  4  grids.  All  of  the  solutions  were  run 
with  the  single  processor,  explicit  time  integration  version  of  Cobalt.  Prior  to  this  study,  we  have 
had  very  little  experience  running  Cobalt  on  which  to  base  a  determination  of  adequate  surface 
grid  resolution.  An  inviscid  study  was  therefore  performed  to  identify  an  adequate  surface  grid 
resolution  to  model  the  important  physics  of  the  problem.  We  encountered  several  difficulties 
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while  generating  viscous  grids  on  the  ONERA  M6  wing  in  TETMESH.  As  a  result,  only  one 
viscous  TETMESH  grid  was  generated.  These  difficulties  are  outlined  below  along  with  a 
presentation  of  the  Cobalt  results. 

4.2.1  Surface  Grid  Resolution  Study 

A  grid  refinement  study  was  performed  to  identify  surface  grid  requirements  for  the 
Cobalt  flow  solver.  Two  grids  (tel  and  te2)  were  generated  with  TETMESH.  A  third  mesh, 
nel,  which  came  from  a  NASA  Langley  workshop  on  unstructured  grid  methods  was  also  used. 
The  grid  sizes  are  summarized  in  Table  4.  The  surface  mesh  and  a  mid-span  cut  through  the 
volume  grid  of  each  mesh  is  shown  in  Figure  65.  Grid  nel  has  a  higher  concentration  of  surface 
grid  nodes  near  the  leading  edge  relative  to  the  other  grids.  The  volume  cell  size  of  the 
TETMESH  grids  grow  very  rapidly  away  from  the  wing  surface.  More  control  over  this 
stretching  would  be  desirable.  Cobalt  Euler  solutions  were  computed  at  run  conditions  of  Mach 
0.84  and  angle  of  attack  of  3.04  degrees.  Predicted  surface  pressures  at  four  spanwise  cuts  are 
compared  with  experimental  data  in  Figtare  66.  The  solutions  on  all  three  meshes  are  similar. 

The  shock  predictions  on  the  te2  mesh  are  sli^tly  crisper  because  of  the  increased  surface 
resolution.  The  greatest  differences  can  be  seen  at  the  65%  span  location  where  the  tel  mesh 
solution  completely  smears  over  the  first  shock.  The  surface  solution  on  mesh  nel  is  smoother 
than  the  solutions  on  meshes  tel  and  te2.  This  may  be  due  to  the  smooth  ordering  of  sairfaces 
faces  in  mesh  nel.  It  appears  that  a  surface  mesh  resolution  somewhere  between  tel  and  te2  is 
required  to  adequately  resolve  the  flowfield  features.  Unfortunately  we  were  unable  to  achieve 
this  level  of  surfece  grid  resolution  in  the  viscous  meshes  we  generated  with  TETMESH.  The 
size  of  our  viscous  TETMESH  grids  was  limited  by  the  memory  requirements  of  TETMESH  and 
the  available  memory  on  the  machine  we  were  running  on. 
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a)  grid  tel  (11,568  surface  faces,  163,930  tetrahedra) 


b)  grid  te2  (32,510  surface  faces,  220,251  tetrahedra) 


c)  grid  nel  (15,279  surface  faces,  231,507  tetrahedra) 


Figure  65.  Surface  Grid  Resolution  Study  -  Tetrahedral  Grids  About  ONERA  M6  Wing. 


Figure  66  Predicted  Surface  Pressure  Coefficient  From  Cobalt  Surface  Grid  Resolution 
Study  on  ONERA  M6  Wing  (Mach  0.84,  a=3.04  degrees). 


4.2.2  Viscous  Results 

TETMESH  was  used  to  generate  prismatic/tetrahedral  grids  for  viscous  analysis  of  the 
ONERA  M6  Wing.  We  were  unable  to  generate  a  viscous  packed  grid  about  an  object  with  a 
symmetry  plane  in  TETMESH.  The  wing  surface  was  therefore  reflected  about  the  symmetry 
plane  and  the  volume  grid  was  generated  on  the  full  wing  geometry  (no  symmetry  plane).  This 
reduced  the  grid  resolution  by  a  factor  of  two  since  our  grid  sizes  were  limited  by  the  amount 
of  available  memory.  The  surface  grid  and  a  cut  through  the  volume  grid  of  a  TETMESH 
generated  viscous  mesh  (tvl)  are  shown  in  Figure  67.  This  grid  contains  20,160  surface  faces 
and  422,614  grid  cells.  The  normal  spacing  at  the  wall  was  selected  to  give  a  y^  value  of  4  at 
a  Re)molds  number  of  1 1.7  million  based  on  the  mean  aerodynamic  chord  length.  A  geometric 
stretching  was  used  to  extrude  the  surface  mesh  across  a  viscous  layer  thickness  which 
corresponds  to  a  y"^  value  of  800.  As  seen  in  Figure  67,  there  is  a  large  discontinuity  in  grid 
spacing  between  the  edge  of  the  extruded  and  tetrahedral  grids.  Extening  the  extruded  region 
further  into  the  field  would  reduce  the  jump  in  grid  spacing,  unfortunately,  attempts  to  generate 
a  grid  with  a  thicker  viscous  layer  resulted  in  crossed  faces  in  the  volume  grid.  Several 
unsuccessful  attempts  were  made  to  eliminate  the  crossed  faces  by  re-meshing  the  surface  grid, 
and  smoothing  the  surface  normals. 

Attempts  were  made  to  run  Cobalt  on  grid  tvl  at  run  conditions  of  Mach  0.7, 0.04  degrees 
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sngle  of  attack,  and  Mach  0.84,  5.06  degrees  angle  of  attack  and  Reynolds  number  of  11.7 
million.  The  Spalart  turbulence  model  was  used  in  all  solutions.  The  solution  at  0.04  degrees 
angle  of  attack  was  run  first-order  for  4,000  iterations  and  second-order  for  22,000  iterations  at 
a  CFL  number  of  0.8.  The  convergence  history  is  shown  in  Figure  68.  The  drag  converges  for 
the  first  10,000  iterations  and  then  begins  to  oscillate.  The  cause  for  this  lack  of  convergence 
is  not  known  but  probably  is  a  consequence  of  the  poor  grid  quality.  The  solution  required  58.7 
hours  of  CPU  time  and  60.9  million  words  of  memory  on  a  CRAY  C90  processor.  Spanwise 
cuts  of  the  surface  pressure  from  this  solution  are  compared  with  experimental  data  in  Figure  69. 
The  comparison  is  relatively  good  with  the  exception  of  the  leading  edge  region  of  the  outboard 
cuts.  The  poor  comparison  in  this  region  may  be  due  to  several  factors  including,  lack  of 
convergence,  lack  of  grid  resolution  in  the  outer  portion  of  the  boundary  layer  and/or  poor  grid 
quality.  The  integrated  forces  from  this  solution  are  shown  in  Table  5.  We  were  unable  to  get 
Cobalt  to  converge  at  5.06  degrees  angle  of  attack.  Attempts  were  made  to  run  Cobalt  on 
grid  wl,  a  tetrahedral  viscous  packed  grid  created  by  the  VGRDD  grid  generator.  This  solution 
would  have  provided  equitable  comparisons  between  the  Cobalt  and  USM3D  flow  solvers  on  the 
same  grid.  Unfortunately,  we  have  been  unable  to  get  a  Cobalt  solution  to  converge  on  this  grid. 


Figure  67.  Surface  and  cut  through  volume  of  viscous  packed  TETMESH  generated  grid 
about  ONERA  M6  wing  (grid  tvl). 
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Figure  68.  Cobalt  solution  convergence  history  for  ONERA  M6  wing  (Mach  0.7,  a=0.04 
degrees,  Re^  =  1 1.7  million,  grid  tvl). 
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Figiire  69.  Cobalt  predicted  surface  pressure  coefficient  at  4  spanwise  cuts  of  ONERA  M6 
wing  (Mach  0.7,  a=0.04  degrees,  Re^  =11.7  million,  grid  tvl). 


4.3  NASTD/MACGS  Results 


A  total  of  3  NASTD  viscous  solutions  were  computed  on  1  MACGS  generated  grid 
(mvl).  For  comparison  purposes,  an  Euler  solution  was  also  computed  on  grid  nel. 

A  multi-zone  hybrid  viscous  packed  grid  was  generated  in  MACGS.  The  surface  and 
symmetry  plane  resolution  are  shown  in  Figure  70.  This  grid  contains  3264  surface  faces, 
73,008  hexahedral  volume  cells,  293,265  tetrahedral  volume  cells,  and  138,148  grid  nodes.  The 
interface  between  the  structured  and  unstructured  grids  is  non-point  matched.  The  resolution  of 
the  unstructured  surface  mesh  at  the  zonal  interface  was  selected  to  minimize  the  jump  in  grid 
spacing  between  the  structured  and  unstructured  meshes.  The  grid  spacing  normal  to  the  surface 
was  selected  to  give  a  y"^  of  4  at  a  Reynolds  number  of  1 1 .7  million.  Care  was  taken  to  place 
the  zone  boundary  between  the  structured  hexahedral  grid  and  the  unstructured  tetrahedral  grid 
far  enough  away  from  the  geometry  to  be  outside  of  the  viscous  dominated  region  of  the  flow. 

NASTD  solutions  were  generated  on  grid  mvl  at  flow  conditions  of  Mach  0.84,  0.04 
degrees  angle  of  attack,  Reynolds  number  of  1 1.7  million,  and  Mach  0.7,  5.06  degrees  angle  of 
attack,  and  Reynolds  number  of  11.7  million.  Turbulence  was  modeled  using  the 
Spalart/Allmaras  (SA)  turbulence  model. 

The  convergence  history  of  the  0.04  degree  angle  of  attack  case  is  shown  in  Figure  71. 
The  solution  was  run  for  12,000  iterations  at  a  CFL  number  of  1.2  although  convergence  was 
achieved  in  roughly  6,000  iterations.  The  solution  required  169.9  hours  of  CPU  time  for  12,000 
iterations  on  an  HP735  workstation  (equivalent  to  48.5  CRAY  C-90  hours)  and  4.9  million  words 
of  memory.  Spanwise  cuts  of  the  predicted  wing  surface  pressure  for  the  0.04  degree  angle  of 
attack  case  are  compared  with  experimental  data  in  Figure  72.  The  surface  pressures  agree  well 
with  the  experimental  data  with  the  exception  of  the  lower  surface  leading  edge  in  the  65%  span 
cut.  The  reason  for  this  discrepancy  is  not  known.  The  integrated  lift  and  drag  from  this 
solution  are  presented  in  Table  5. 

Surface  flow  patterns  from  the  NASTD  solution  for  the  5.06  degree  angle  of  attack  case 
are  shown  in  Figure  73.  These  patterns  reveal  a  substantial  shock-induced  separation  on  the 
outboard  portion  of  the  wing.  Spanwise  cuts  of  the  wing  surface  pressure  are  compared  with 
experimental  data  in  Figure  74.  The  surface  pressures  agree  very  well  with  the  experimental 
data.  The  predicted  shock  location  matches  the  experimental  data  indicating  that  the  location  and 
of  the  separation  is  well  predicted.  There  is  some  discrepancy  between  the  NASTD 
solution  and  the  experimental  data  just  aft  of  the  separation  line.  This  may  be  due  to  the 
turbulence  model. 

To  investigate  the  effects  of  the  turbulence  model  choice,  the  5.06  degree  angle  of  attack 
case  was  recomputed  in  NASTD  using  the  SST  turbulence  model  (5).  This  turbul^ce  model  is 
a  two-equation  k-co  turbulence  model.  The  convergence  of  this  solution  was  similar  to  the 
Spalart  Allmaras  turbulence  model  solution,  however,  the  solution  time  per  iteration  was  lower 
due  to  fewer  sub-iterations  required  by  the  turbulence  model.  The  predicted  surface  pressure 
from  this  solution  is  compared  with  the  Spalart  turbulence  model  solution  and  experimental  data 
in  Figure  75.  The  two-equation  turbulence  model  appears  to  have  only  minor  effects  on  the 
solution.  Most  of  the  differences  are  in  the  separation  region  where  the  two-equation  model 
solution  is  marginally  better. 
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Pressure-Force  Coefficient 


Figure  72.  NASTD  predicted  surface  pressure  coefficient  at  4  spanwise  cuts  of  ONERA  M6 
wing  (Mach  0.7,  a=0.04  degrees,  Re^  =11.7  million,  grid  mvl). 


Figure  73.  NASTD  predicted  particle  traces  confined  to  surface  of  ONERA  M6  wing  (Mach 
0.84,  a=5.06  degrees,  Reg=11.7  million,  grid  mvl). 
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Figure  74.  NASTD  predicted  surfece  pressure  coefficient  at  6  spanwise  cuts  of  ONERA  M6 
wing  (Mach  0.84,  a=5.06  degrees,  Re^  =11.7  million,  grid  mvl). 
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Figure  75.  Effect  of  turbulence  model  on  NASTD  predicted  surface  pressure  coefficient  for 
ONERA  M6  wing  (Mach  0.84,  a=5.06  degrees,  Re,  =  11.7  million,  grid  mvl). 


4.4  USM3DA^GRID  Results 

A  total  of  3  USM3D  viscous  solutions  were  computed  on  2  VGRID  generated  grids.  For 
comparison  purposes,  an  Euler  solution  was  also  computed  on  grid  nel.  The  viscous  USM3D 
computations  were  performed  on  two  existing  grids,  wl  and  vv2.  A  summary  of  the  grid 
charactenstics  for  each  grid  is  presented  in  Table  4.  The  surface  grid  and  symmetry  plane  of 
gnds  wl  and  w2  are  shown  in  Figure  76.  The  surface  mesh  is  the  same  for  both  grids  and  is 


49 


stretched  anisotropically  in  the  streamwise  direction  in  order  to  reduce  the  overall  grid  size.  The 
viscous  layer  merges  smoothly  into  the  inviscid  grid  with  a  minimal  jump  in  grid  spacing.  The 
grid  spacing  of  the  first  point  normal  to  the  wall  results  in  a  y"^  value  of  80  for  grid  wl  and  4 
for  grid  vv2  at  a  Reynolds  number  of  1 1 .7  million.  The  large  y"^  value  of  grid  wl  necessitates 
the  use  of  a  wall  fimction  to  model  the  near  wall  viscous  effects  while  grid  w2  was  used  without 
a  wall  function. 

Two  solutions  were  run  on  grid  wl  using  the  Spalart/Allmaras  turbulence  model  and  a 
wall  function.  The  solutions  were  run  at  flow  conditions  of  Mach  0.847,  Re=l  1.7  ntillion  and 
angles  of  attack  of  0.04  and  5.06  degrees.  The  0.04  degree  angle  of  attack  case  was  run  for  3600 
iterations  at  a  CFL  number  of  200.  The  solution  required  29.4  hours  of  CPU  time  and  63.4 
million  words  of  memory  on  a  CRAY  C-90  processor.  Spanwise  cuts  of  the  predicted  wing 
surface  pressure  for  the  0.04  degree  angle  of  attack  case  are  compared  with  e3q)erimental  data 
in  Figure  77.  The  surface  pressures  agree  well  with  the  experimental  data.  The  integrated  lift 
and  drag  firom  this  solution  are  presented  in  Table  5. 

Surface  flow  patterns  fijom  the  USM3D  solution  for  the  5.06  degree  angle  of  attack  case 
are  shown  in  Figure  78.  As  with  the  NASTD  solution,  a  large  separation  region  is  predicted  on 
the  outboard  portion  of  the  wing,  however,  the  size  and  location  of  this  separation  are  different. 
Spanwise  cuts  of  the  wing  surface  pressure  are  compared  with  experimental  data  in  Figure  79. 
The  solution  at  cuts  inboard  of  the  separation  agree  well  with  the  experimental  data,  however, 
the  solution  at  the  cuts  outboard  of  the  separation  show  that  the  separation  location  is  predicted 
too  far  forward.  This  may  be  due  to  the  use  of  a  wall  function  at  a  separated  flow  condition. 

The  effect  of  the  wall  function  on  the  solution  was  investigated  by  computing  a  full 
viscous  (no  wall  function)  USM3D  solution  on  grid  vv2.  The  surface  flow  patterns  fi:om  this 
solution  are  shown  in  Figure  80  and  the  predicted  surface  pressure  is  shown  in  Figure  81.  Once 
again  the  surface  flow  traces  indicate  a  substantial  shock  induced  separation  on  the  outboard 
portion  of  the  wing.  Without  the  wall  function,  the  surface  pressure  agrees  much  closer  with  the 
etqjerimental  data  at  the  outboard  cut  locations. 
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Figure  76.  Mesh  on  surface  and  symmetry  plme  of  two  VGRID  generated  tetrahedral  grids 
about  the  ONERA  M6  wing  (grids  wl  and  w2). 


wing  (Mach  0.7,  a=0.04  degrees,  Re^  =11.7  million,  grid  wl). 
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Figure  78.  USM3D  predicted  particle  traces  confined  to  surface  of  ONERA  M6  wing  (Mach 
0.84,  a=5.06  degrees,  Rej=11.7  million,  grid  wl). 
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Figure  79.  USM3D  predicted  surface  pressure  coefficient  at  6  spanwise  cuts  of  ONERA  M6 
wing  (Mach  0.84,  a=5.06  degrees,  Re,.  =  11.7  million,  grid  wl). 
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Figure  80.  USM3D  predicted  particle  traces  coniined  to  surface  of  ONERA  M6  wing  (Mach 
0.84,  a=5.06  degrees,  Rec=11.7  million,  grid  w2). 
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Figure  81.  USM3D  predicted  surface  pressure  coefficient  at  6  spanwise  cuts  of  ONERA  M6 
wing  (Mach  0.84,  a=5.06  degrees,  Re^  =11.7  million,  grid  w2). 


4.5  Flow  Solver  Comparisons 

In  this  section  the  efficiency  and  accuracy  of  the  ONERA  wing  flow  solutions  from  the 
three  methods  are  compared.  A  summary  of  all  ONERA  wing  flow  solutions  generated  during 
this  study  is  given  in  Table  5. 
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Table  5.  ONERA  Wing  Run  Matrix 
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The  different  grid  topologies  used  with  each  flow  solver  for  viscous  computations 
precludes  an  equitable  one-to-one  comparison  of  flow  solver  accuracy  and  efficiency  independent 
of  the  grid.  However,  such  a  comparison  can  be  made  for  inviscid  flow  since  each  of  the  three 
methods  can  compute  Euler  solutions  on  tetrahedral  grids.  For  this  purpose,  Euler  solutions  were 
computed  with  all  three  methods  on  grid  nel  at  run  conditions  of  Mach  0.84,  and  3.06  degrees 
angle  of  attack.  In  Figure  82,  the  predicted  surface  pressure  from  the  three  methods  is  compared 
with  experimental  data.  All  three  solutions  agree  relatively  well  with  the  experimental  data.  The 
Cobalt  solution  appears  to  be  marginally  better  than  the  other  two  methods  particularly  near  the 
suction  peaks.  However,  the  Cobalt  solution  also  had  the  longest  computer  run  time  of  the  three 
methods.  The  memory  and  CPU  requirements  for  each  of  the  three  methods  are  tabulated  in 
Table  5.  Both  USM3D  and  Cobalt  are  cell  centered  methods  and  require  approximately  600% 
more  memory  than  the  nodal  based  algorithm  used  in  NASTD.  The  CPU  requirements  of  Cobalt 
were  significantly  larger  than  USM3D  and  NASTD  primarily  due  to  the  large  number  of 
iterations  required  to  reach  convergence.  NASTD  required  the  least  amount  of  CPU  of  the  three 
methods  most  likely  due  to  the  lower  amount  of  work  required  by  the  nodal  based  scheme. 

The  results,  obtained  on  the  ONERA  wing  at  Mach  0.84,  Re=11.7  million,  and  0.04 
degrees  angle  of  attack,  provide  the  opportunity  to  compare  the  viscus  prediction  capability  of 
the  flow  solvers  in  the  absence  of  any  flow  separation.  The  predicted  surface  pressures  from 
each  of  the  methods  is  compared  with  the  experimental  data  in  Figure  83.  Without  any  flow 
separation  each  of  the  three  methods  does  a  fairly  good  job  of  predicting  the  surface  pressures. 
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TTie  Cobalt  solution  is  the  worst  of  the  three  probably  due  to  the  poor  grid  qxiality.  The 
mtegrated  drag  from  each  of  the  three  methods  is  shown  in  Table  5.  The  predicted  pressure  drag 
from  the  3  methods  ranges  from  0.0018  to  0.0103.  Likewise  the  skin  friction  drag  varies  from 
.0035  to  0.0296.  Unfortunately,  without  running  the  flow  solvers  on  the  same  grid  it  is 
impossible  to  determine  if  this  discrepancy  is  due  to  grid  differences  or  to  differences  in  the  flow 
solver  algorithms.  The  memory  requirements  of  USM3D  and  Cobalt  were  similar  at  60.9  MW 
for  Cobalt  and  63.4  MW  for  USM3D.  The  NASTD  memory  requirements  were  much  lower  at 
4.9  MW.  This  is  due  to  several  factors  including  more  efficient  storage  of  structured  grid  data, 
nodal  based  scheme,  and  multi-zone.  The  run  times  for  NASTD  (after  6,000  iterations)  and 
USM3D  were  roughly  equivalent  at  103.6  HP  CPU  hours  (equivalent  to  29.6  CRAY  C-90  CPU 
hours)  and  29.39  hours  of  CRAY  C-90  CPU  time.  The  run  time  for  Cobalt  was  almost  twice 
as  long  at  58.67  CRAY  CPU  hours. 

The  second  set  of  run  conditions,  Mach  0.84,  Re=11.7  million,  5.06  degrees  angle  of 
attack,  provides  a  comparison  of  the  different  methods  in  a  turbulent  separated  flow.  In  Figure 
84,  spanwise  cuts  of  the  predicted  surface  pressure  from  the  NASTD  and  USM3D  (without  wall 
ftmction)  solutions  are  compared  with  experimental  data.  Considering  the  complexity  of  the 
flowfield,  both  solutions  agree  remarkably  well  with  the  experimental  data.  The  NASTD  results 
predict  the  shock  boundary  layer  interaction  slightly  better  showing  slightly  more  smearing  across 
the  shock  than  the  USM3D  results.  The  NASTD  results  also  match  the  experimental  data  slightly 
better  in  the  separated  regions  of  the  flowfield.  These  differences  could  be  due  to  the  different 
grid  resolutions  of  the  two  solutions.  USM3D  required  98.2  MW  of  memory  compared  to  5.7 
MW  for  NASTD.  The  run  times  of  the  two  methods  were  46.6  CRAY  C-90  CPU  hours  for 
USM3D  and  220.99  HP  735  CPU  hours  (equivalent  to  63.14  CRAY  C-90  CPU  hours)  for 
NASTD. 
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Figure  82.  Comparison  of  Euler  flow  solutions  on  ONERA  M6  wing  (Mach  0.84,  a=3.04 
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Figure  84.  Comparison  of  Navier-Stokes  solutions  on  ONERA  M6  wing  (Mach  0.84,  a=5.06 

degrees,  Ree=11.7  million) 
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5.  METHOD  ASSESSMENT 


5.1  Accuracy 

A  direct  comparison  of  the  relative  accuracy  of  the  three  methods  evaluated  in  this  study 
is  not  possible  since  each  method  utilizes  a  different  grid  cell  topology.  The  one  direct 
comparison  that  was  made,  an  Euler  solution  on  the  ONERA  wing,  did  not  show  sigmficant 
differences  in  the  accuracy  of  the  three  methods.  There  was  some  variation  in  the  acctiracy  of 
the  three  methods  for  viscous  computations.  The  accuracy  of  Cobalt  and  NASTD  was  very 
similar  for  flat  plates  with  Cobalt  having  slightly  higher  accuracy.  The  USM3D  solutions  without 
the  wall  function  were  the  worst  of  the  three  methods.  However,  the  USM3D  solutions  with  a 
wall  function  offered  the  best  accuracy.  These  results  are  inconsistent  with  previous  results 
generated  by  the  code  author  (1).  Further  study  is  required  to  determine  the  cause  of  this 
discrepancy. 

All  of  the  viscous  unstructured  grid  methods  evaluated  in  this  report  have  the  ability  to 
provide  accuracy  approaching  that  of  existing  structured  grid  CFD  methods  provided  adequate 
grid  resolution  and  quality  is  employed.  However,  unstructured  grid  resolution  and  quality 
requirements  are  less  well  known  than  for  structured  grids  adding  an  element  of  risk  to  the 
solution.  The  grid  requirement  guidelines  developed  in  this  report  provides  a  starting  point  for 
the  grid  developer.  Additional  grid  requirements  for  other  gnd  types  and  mixtures  of  grid  types 
need  to  be  developed.  In  addition  the  effects  of  grid  cell  shape  on  solution  accuracy  for  each 
flow  solver  need  to  be  identified.  Once  these  attributes  are  known,  mixed  cell  grids  such  as 
those  utilized  by  TETMESH/Cobalt  should  provide  much  better  control  over  grid  quality  yielding 

a  distinct  advantage  over  the  tetrahedral  type  grids. 

The  aspect  of  grid  qtiality  was  not  explored  in  this  study  but  can  have  a  substantial  effect 
on  the  solution  accuracy.  Guidelines  need  to  be  developed  to  aid  the  end  user  in  ass^sing  the 
quahty  of  the  grid.  This  is  particularly  important  for  unstructured  grids  since  it  is  almost 
impossible  to  asses  the  grid  quality  visually.  These  guidelines  should  eventually  be  used  in 
concert  with  an  unstructured  grid  adaptation  method  to  remove  the  burden  of  grid  quality 
assessment  from  the  user  by  automatically  assuring  the  grid  meets  the  quality  and  resolution 

guidelines. 


5.2  Efficiency 

A  large  variation  in  flow  solver  efficiency  was  identified  in  this  study.  This  was  primarily 
a  consequence  of  the  different  time  integration  and  discretization  schemes  used  by  each  metl^. 
Implicit  time  integration  was  found  to  be  necessary  for  adequate  solution  convergence.  The 
exphcit  time  integration  scheme  in  Cobalt  was  found  to  be  unacceptably  slow  for  viscous 
dominated  flows  since  the  low  CFL  numbers  severely  limited  the  convergence  rate  m  the 
boundary  layer.  The  implicit  schemes  used  in  NASTD  and  USM3D  provided  much  better 
convergence.  The  point  Jacobi  implicit  scheme  in  USM3D  provided  the  fastest  convergence  of 

the  three  methods  at  the  expense  of  more  memory  usage. 

The  memory  requirements  of  Cobalt  and  USM3D  were  both  very  large  primarily  due  to 
their  unstructured  data  structures  and  the  use  of  cell  based  schemes.  USM3D  has  the  Imgest 
memory  requirements  of  the  three  methods.  NASTD  utilized  memory  much  more  efficiently 
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through  the  use  of  structured  data  structures  in  the  viscous  layer  and  a  nodal  based  scheme. 
Future  use  of  these  methods  on  complex  geometries  may  be  constrained  by  their  memory  usage. 
The  Cobalt  and  USM3D  ONERA  wing  solutions  run  in  this  study  approached  the  memory  limits 
of  currently  available  machines. 

The  largest  contributor  to  the  cycle  time  of  a  CFD  solution  is  in  the  grid  development 
time.  The  unstructured  grids  used  by  these  methods  offer  the  potential  to  drastically  reduce  this 
time  over  existing  structured  grid  methods.  However,  current  unstructured  grid  generation 
systems  will  require  a  great  deal  of  development  before  they  can  match  the  capability  and  ease 
of  use  of  existing  structured  grid  generation  systems.  MACGS  provided  the  largest  array  of 
capabilities  and  ease  of  use  of  the  grid  generation  tools  evaluated  in  this  study.  Unfortunately 
this  system  relies  on  zonal  hybrid  grid  technology.  This  technology  reduces  the  efficiency  of  the 
grid  generation  process  since  separate  structured  and  unstructured  grid  blocks  must  be  developed 
around  the  geometry.  The  VGRID  tool  provided  the  fastest  griding  times  from  start  to  finish  of 
the  evaltiated  packages.  The  mixed  cell  grids  used  by  TETMESH  should  provide  the  greatest 
flexibility  to  the  grid  developer  and  thereby  the  fastest  grid  generation  times.  Unfortunately  this 
efficiency  was  not  realized  because  of  usability  issues  discussed  below. 

53  Ease  of  Use 

One  of  the  most  important  attributes  of  a  CFD  solution  method  is  its  ease  of  use.  Ease 
of  use  is  a  subjective  issue.  The  primary  virtues  that  contributed  to  ease  of  use  in  this  study 
included  quick  learning  curve,  documentation,  robustness,  abundance  of  useful  features,  and 
intuitiveness  of  the  approach.  We  had  a  great  deal  of  expertise  with  the  MACGS/NASTD 
software  package  prior  to  this  study.  This  biases  our  ability  to  make  an  impartial  assessment  of 
the  ease  of  use  of  these  packages.  MACGS/NASTD  are  therefore  not  evaluated  for  ease  of  use 
in  this  study. 

Surprisingly  few  difficulties  were  encountered  during  this  study  given  the  relative 
immaturity  of  the  methods.  The  few  problems  that  were  encoimtered  usually  happened  when 
working  outside  the  environment  for  which  the  method  had  been  previously  tested  by  the 
developers.  Two  examples  are  running  USM3D  widiout  a  wall  function  and  running  Cobalt  on 
tetrahedral  viscous  grids.  Continued  development  of  these  solvers  should  improve  their 
applicability  to  a  wider  range  of  problems.  Leammg  to  use  each  of  the  flow  solvers  was  very 
straightforward.  The  Cobalt  docrunentation  was  excellent  and  provided  a  complete  description 
of  all  features.  USM3D  did  not  come  with  any  documentation  other  than  the  comments  in  the 
Fortran  source  code.  A  users  manual  would  be  a  more  convenient  aid  to  the  new  user.  Both 
Cobalt  and  USM3D  are  relatively  new  software  tools  and  therefore  do  not  have  the  large  number 
of  features  and  capabilities  such  as  multiple  time  and  spatial  niimerical  schemes,  large  selection 
of  boundary  conditions,  dissipation  operators,  or  parallel  processing  cjpability  that  can  be  foimd 
on  more  mature  CFD  solvers.  As  development  of  these  packages  continues  the  number  of 
features  will  continue  to  grow. 

The  greatest  obstacle  to  a  new  user  of  any  of  these  methods  is  the  large  learning  curve 
for  the  grid  generation  tools.  All  of  the  grid  generation  packages  required  several  weeks  of  study 
to  achieve  a  basic  level  of  knowledge.  Documentation  when  available  aided  in  the  learning 
process.  However,  all  of  the  grid  generation  tools  could  use  improvements  to  the  documentation. 
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The  GRIDTOOLA^ GRID  documentation  provided  a  good  description  of  GRIDTOOL  but  did  not 
include  any  mformation  on  VGRED  or  POSTGRID.  The  TETMESH  documentation  was  very 
sparse.  A  complete  description  of  all  of  the  features  in  TETMESH  with  examples  of  their  use 
would  aid  the  new  user  considerably. 

Both  the  GREDTOOLA^GRID  and  TETMESH  packages  are  relatively  immature  and  as 
such  do  not  include  a  large  array  of  features  that  are  essential  in  the  production  CFD 
environment.  Capabilities  such  as  direct  links  to  CAD  packages,  ability  to  modify/repair 
geometry,  complete  control  of  grid  distributions  on  the  surface  and  in  space  were  noticeably 
missing.  Both  of  these  packages  also  will  only  run  on  Silicon  Graphics  workstations  which  limits 
their  usefulness. 

Of  the  two  packages  the  GRIDTOOLA^GRID  offered  the  most  capability.  Once  the 
process  was  learned,  the  sequence  of  events  to  generate  a  grid  in  GRIDTOOL A^GRID  was 
intuitive  and  relatively  strai^t  forward  with  the  exception  of  the  jumps  required  between  the 
GRIDTOOL,  VGRID,  and  POSTGRID  software  tools.  The  user  has  good  control  over  the  grid 
distribution  through  the  specification  of  point  and  line  sources.  The  quality  of  the  resulting  grids 
was  found  to  be  very  good.  A  unique  restart  capability  in  VGRID  allows  very  large  grids  to  be 
developed  on  relatively  small  memory  machines. 

The  TETMESH  software  uses  a  unique  approach  to  generate  grids  containing  a  mixture 
of  prismatic  and  tetrahedral  cells.  The  TETMESH  grid  generation  process  is  straight  forward  for 
tetrahedral  grids.  For  extruded  prismatic  grids  however,  the  process  was  not  intuitive  requiring 
a  number  of  steps  to  be  performed  in  a  specific  order.  Without  documentation  or  online  help  it 
was  very  difficult  to  determine  the  ^propriate  steps. 

During  the  cotu^e  of  this  study  several  limitations  were  encountered  that  made  TETMESH 
very  difficult  to  use.  The  user  does  not  have  adequate  control  over  the  grid  distributions  on  the 
surface  mesh  and  almost  no  control  over  the  distribution  of  the  volume  mesh.  Boundary 
conditions  cannot  be  specified  or  changed  on  the  surfaces  once  the  grid  has  been  generated. 
TETMESH  was  fotmd  to  require  a  large  amount  of  memory  and  there  were  several  instances 
when  the  program  bombed  while  generating  a  mesh.  The  surface  extrusion  capability  does  not 
work  in  conjunction  with  a  symmetry  plane  requiring  the  entire  geometry  to  be  modeled  thereby 
exacerbating  the  large  memory  usage  problem.  The  surface  extrusion  capability  is  very  fragile. 
Several  items  can  cause  the  extrusion  to  fail  including  a  convex  geometry  or  poor  surface  mesh. 
The  tools  that  are  available  to  help  with  these  problems,  such  as  normal  smoothing,  did  not  seem 
to  help  in  the  problems  we  encountered. 
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6.  CONCLUSIONS/RECOMMENDATIONS 


A  systematic  assessment  of  the  strengths  and  weaknesses  of  three  viscous  flowfield 
simulation  methods  for  unstructured  grids  has  been  performed.  The  methods  evaluated  rq)resent 
the  primary  viscous  unstractured  griding  approaches  currently  in  use  including  tetrahedral,  zonal 
hybrid,  and  grids  containing  a  mixture  of  tetrahedral  and  prismatic  grids.  The  strengths  and 
weaknesses  of  each  method  have  been  identified  through  a  systematic  study  of  the  flowfields 
about  a  flat  plate  and  ONERA  M6  wing.  Solutions  were  generated  over  a  range  of  run 
conditions  allowing  comparisons  of  a  wide  variety  of  flow  field  features  including  transonic  and 
separated  flows.  Each  of  the  methods  were  evaluated  based  on  their  relative  efficiency,  accuracy, 
ease  of  use,  and  robustness.  Predictions  of  skin  friction,  velocity  profiles,  and  surface  pressure 
with  theoretical  and  experimental  data  were  used  to  evaluate  solution  accuracy.  Flow  solver 
efficiency  was  evaluated  by  comparing  memory  and  processor  requirements  of  each  method. 

The  different  approaches  used  by  each  of  the  software  packages  resulted  in  a  unique  set 
of  advantages  over  the  other  methods.  The  mixed  cell  unstmctured  grid  used  by 
Cobalt/TETMESH  provides  a  great  deal  of  flexibility  which  has  the  potential  to  reduce  grid 
generation  cycle  time  and  improve  solution  quality  through  improved  grid  quality.  The 
MACGS/NASTD  software  package  was  found  to  be  the  most  efficient  of  the  methods.  The  use 
of  structured  grid  data  structures  and  a  nodal  based  scheme  gave  NASTD  the  lowest  memory 
usage  of  the  methods  and  solution  run  times  about  equal  to  USM3D  and  almost  half  that  of 
Cobalt.  USM3D  utilized  a  point  Jacobi  integration  scheme  in  conjimction  with  wall  functions 
to  provide  the  fastest  run  times  of  the  three  methods.  The  GRIDTOOLA^ GRID  grid  generation 
package  was  the  easiest  to  use  of  the  grid  generation  systems  and  provided  the  fastest  grid 
generation  cycle  time. 

Each  of  the  methods  were  foxmd  to  provide  adequate  accuracy  for  engineering  analysis 
provided  adequate  grid  resolution  and  quality  are  employed.  A  set  of  grid  resolution  guidelines 
was  developed  for  grid  wall  spacing  and  number  of  points  in  the  boundary  layer.  Additional 
guidelines  are  needed  to  define  required  grid  resolution  and  quality  for  specific  types  of  flow 
fields.  In  addition  grid  adaptation  techniques  need  to  be  developed  to  automatically  ensure  grid 
quality  and  resolution. 

Several  improvements  can  be  made  to  improve  the  usability  of  these  unstructured  grid 
methods.  The  biggest  roadblock  to  future  use  of  these  systems  is  in  the  grid  generation  systems. 
The  primary  focus  of  future  development  should  therefore  be  in  improving  these  systems  by 
developing  better  documentation  and  improved  ease  of  use.  Future  flow  solver  method 
developments  should  focus  on  reducing  memory  requirements  and  improving  efficiency. 
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APPENDICES 


Residual 


Cobalt  Flat  Plate  Solutions,  Mach  0.1,  Re=1.77x^0^  Turbulent  Flow 


Figure  A-11.  Cobalt  Residual  History,  Grid  c-6  Figure  A-12.  Cobalt  Force  Convergence  History,  Grid  c-6 


Iterations 


Figure  A-13.  Cobalt  Residual  History,  Grid  c-7  Figure  A-14.  Cobalt  Force  Convergence  History,  Grid  c-7 


Figure  A-15.  Cobalt  Residual  History,  Grid  c-8  Figure  A-16.  Cobalt  Force  Convergence  History,  Grid  c-8 
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Figure  A-21.  Cobalt  Residual  History,  Grid  c-14  Figure  A-22.  Cobalt  Force  Convergence  History,  Grid  c-14 


Cobalt  Flat  Plate  Solutions,  Mach  0.7,  Re=12xl0^  Turbulent  Flow 
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Figure  A-29.  Cobalt  Residual  History,  Grid  c-9  Figure  A-30.  Cobalt  Force  Convergence  History,  Grid  c-9 


Figure  A-31.  Cobalt  Residual  History,  Grid  c-1 1  Figure  A-32.  Cobalt  Force  Convergence  History,  Grid  c-1 1 


Appendix  B:  NASTD  Solution  Convei^ence  Plots 
NASTD  Flat  Plate  Solutions,  Mach  0.1,  Re=1.77xl0‘,  Laminar  Flow 
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Figure  B-3.  NASTD  Residual  History,  Grid  n-9  Figure  B-4.  NASTD  Force  Convergence  History,  Grid  n-9 
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Figure  B-5.  NASTD  Residual  History,  Grid  n-1 1  Figure  B-6.  NASTD  Force  Convergence  History,  Grid  n-1 1 
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Figure  B-7.  NASTD  Residual  History,  Grid  n-12  Figure  B-8.  NASID  Force  Convergence  History,  Grid  n-12 


Figure  B-9.  NASTD  Residual  History,  Grid  n-13  Figure  B-10.  NASTD  Force  Convergence  History,  Grid  n-13 
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Figure  B-13.  NASTD  Residual  History,  Grid  ii-6  Figure  B-14.  NASTD  Force  Convergence  History,  Grid  n-6 


Figure  B-15.  NASTD  Residual  History,  Grid  n-7  Figure  B-16.  NASTD  Force  Convergence  History,  Grid  n-7 


72 


4000 


Iterations 


1COOO 


24000 


Iteraiions 


Figure  B-17.  NASTD  Residual  History,  Grid  n-8  Figure  B-18  JSfASTD  Force  Convergence  History,  Grid  n-8 
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Figure  B-19.  NASTD  Residual  History,  Grid  n-9  Figure  B-20.  NASTD  Force  Conveigence  History,  Grid  n-9 
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Figure  B-21.  NASTD  Residual  History,  Grid  n-11  Figure  B-22.  NASTD  Force  Convergence  History,  Grid  n-1 1 


NASTD  Flat  Plate  Solutions,  Mach  0.7,  Re=12xl0‘,  Turbulent  Flow 


Figure  B-29.  NASTD  Residual  History,  Gnd  n-7  Figure  B-30.  NASTD  Force  CoQveigence  History,  Grid  n-7 
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Appendix  C:  USM3D  Solution  Convei^ence  Plots 
USM3D  Flat  Plate  Solutions,  Mach  0.1,  Re=1.77xl0^  T^aminar  Flow 
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Figure  C-1.  USM3D  Convergence  History,  Grid  u-3  Figure  C-2.  USM3D  Convergence  History,  Grid  u-6 

USM3D  Flat  Plate  Solutions,  Mach  0.1,  Re=1.77xl0‘,  Turbulent  Row 


Figure  C-3.  USM3D  Convergence  History,  Grid  u-2  Figure  C-4.  USM3D  Convergence  History,  Grid  u-3 


Figure  C-5.  USM3D  Convergence  History,  Grid  u-4  Figure  C-6.  USM3D  Convergence  History,  Grid  u-5 
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USM3D  Flat  Plate  Solutions,  Mach  0.7,  Re=12xl0‘,  Turbulent  Flow 
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Figure  C-7.  USM3D  Convergence  History,  Grid  u-1  Figure  C-8.  USM3D  Convergence  History.  Grid  u-3 


Figure  C-9.  USM3D  Ccmvergence  History,  Grid  u-6 
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